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Abstract 

Properties of six- component electromagnetic field solutions of a matrix form of the 
Maxwell equations, analogous to the four-component solutions of the Dirac equation, are 
described. It is shown that the six-component equation, including sources, is invariant un- 
der Lorentz transformations. Complete sets of eigenf unctions of the Hamiltonian for the 
electromagnetic fields, which may be interpreted as photon wave functions, are given both 
for plane waves and for angular-momentum eigenstates. Rotationally invariant projection 
operators are used to identify transverse or longitudinal electric and magnetic fields. For 
plane waves, the velocity transformed transverse wave functions are also transverse, and the 
velocity transformed longitudinal wave functions include both longitudinal and transverse 
components. A suitable sum over these eigenf unctions provides a Green function for the 
matrix Maxwell equation, which can be expressed in the same covariant form as the Green 
function for the Dirac equation. Radiation from a dipole source and from a Dirac atomic 
transition current are calculated to illustrate applications of the Maxwell Green function. 



1. Introduction 

For quantum mechanics to provide a complete description of nature, it is necessary to 
have a wave function for something as important as electromagnetic radiation or photons. 
This has become increasingly relevant as the number of experiments on single photon pro- 
duction and detection, motivated by interest in the fields of quantum computation and 
quantum cryptography, has grown rapidly over the past two decades fl]. The history of 
theoretical efforts to define photon wave functions dates back to the early days of quantum 
mechanics and is still unfolding. Overviews have been given in 0,0, HI- However, there is 
not yet a consensus on the form a photon wave function should take or the properties it 
should have. Further investigation of these questions is warranted, and possible answers are 
given in this paper. 

Quantum electrodynamics (QED) accurately describes the interaction of radiation with 
free electrons and electrons bound in atoms, but as it is formulated in terms of an S matrix, 
asymptotic states, and Feynman diagrams, it does not readily lend itself to the description of 
the time evolution of radiation. In particular, interference effects or the space-time behavior 
of a photon wave packet would be more naturally described in the framework of wave 
mechanics with a wave function for a single photon. 
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There are a number of requirements that need to be be imposed on a formalism for 
a quantum mechanical description of photons. First, the predicted behavior of radiation 
should be consistent with the Maxwell equations. Time-dependent solutions of the Maxwell 
equations provide the basis for both classical electromagnetic theory and QED, and it can 
be expected that a photon wave function should also be based on solutions of the Maxwell 
equations. This means that the wave function is simultaneously the solution of both of the 
first-order Maxwell equations with time derivatives and not just a solution of a second-order 
scalar wave equation. 

A second requirement is that the wave functions obey the quantum mechanical principle 
of linear superposition. Simply stated, this means that if two wave functions describe possible 
states of radiation, then a linear combination of these wave functions also describes a possible 
state. For example, a wave function for circularly polarized radiation can be written as a 
linear combination of two wave functions for linearly polarized radiation, all of which must 
be solutions of the same wave equation. 

Another requirement is that the formalism be Lorentz invariant in order to properly 
describe the space-time behavior of radiation. An approximation scheme like the reduction 
of the Dirac equation to obtain the Schrodinger equation for electron velocities that are 
small compared to the speed of light is not an option for radiative photons. 

Finally, it is necessary for the formalism to provide the tools for methods associated with 
quantum mechanics. This includes a wave equation with a Hamiltonian that describes the 
time development of states, wave functions that comprise a complete set of eigenfunctions of 
the Hamiltonian, normalizable states with a probability distribution that corresponds to the 
location of the photon, a law of conservation of probability, operators with expectation values 
for observables, and wave packets that realistically describe the propagation of photons in 
space and time. 

To arrive at a wave equation that addresses these requirements, we examine an approach 
in which the four-component matrix Dirac equation for a spin one-half electron is adapted 
to a six-component form of the Maxwell equations for a spin-one photon. This version of 
the Maxwell equations is a direct extension of the Dirac equation for the electron in which 
two-by-two Pauli matrices are replaced by analogous three-by-three matrices. Since the 
quantum mechanical properties of the Dirac equation, Hamiltonian, and wave functions are 
well understood and tested experimentally, it is natural to consider the analogous Maxwell 
equation, Hamiltonian, and wave functions as a quantum mechanical description of photons. 

There are fundamental differences between the Dirac equation and the matrix Maxwell 
equation, so the extension requires a detailed analysis. The most prominent difference is the 
fact that there is a possible source term in the Maxwell equation which has no analog for 
the Dirac equation [5( . Also, some properties of the three- by-three spin matrices differ from 
those of the Pauli matrices, even though they have the same commutation relations. 

Linear operators that are representations of the inhomogeneous Lorentz group can replace 
the wave equation of a system for free electrons or transverse photons with no sources @, 0, S] , 
but the source terms and longitudinal solutions of the Maxwell equation fall outside this 
framework. Taking the view that the Maxwell equation with a source is the most direct 
contact with experiment, our approach is to start from the matrix Maxwell equation with 
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a source term and explicitly work out the Lorentz transformations of the solutions. It is 
shown that the six-component equation is invariant under Lorentz transformations, as it 
should be, but this is not self-evident, since the source term is essentially the three- vector 
current density. 

Next, six-component solutions are constructed and shown to be complete sets of orthogo- 
nal coordinate-space eigenf unctions of the Maxwell Hamiltonian, parameterized by physical 
properties, such as linear momentum, angular momentum, and parity. These properties 
are associated with operators that commute with the Hamiltonian. Complete sets of both 
plane-wave solutions and angular-momentum eigenfunctions are given. Bilinear products of 
normalizable linear combinations of these functions provide expressions for the probability 
density and flux. The eigenfunctions are further classified according to whether they rep- 
resent transverse or longitudinal states. These properties are associated with the electric 
and magnetic fields, with the result that under a velocity boost, the transformed transverse 
solutions are also transverse, unlike solutions corresponding to a transverse vector poten- 
tial. Moreover, by summing over both transverse and longitudinal solutions, we obtain a 
covariant Green function for the Maxwell equation, which is of the same form as the Green 
function for the Dirac equation. 

Solutions are obtained directly from the Maxwell equation, with no recourse to a vector 
potential. This avoids problems such as extra polarization components and ambiguities 
associated with gauge transformations j§, [hJ HI- Although integrals over a closed path 
of the potential may be observables, as discussed in [l2j, such integrals can be expressed 
in terms of the magnetic flux through the loop 13j, so it is expected that the fields alone 
provide a complete description of electrodynamics. It is also clear that photon wave functions 
are closely aligned with electric and magnetic fields, and an approach that starts with 
classical electrodynamics expressed in terms of fields only provides a natural framework for 
the transition to the wave mechanics of photons. A possible advantage of using a vector 
potential is that it is the solution of a scalar wave equation, which has a well-known Green 
function. However, this advantage is offset by the fact that we provide a covariant Green 
function for the Maxwell equation. 

This paper is organized as follows. In Sec. [2] the vector Maxwell equations and the Dirac 
equation are stated to define notation. The algebra of three-component spin matrices is 
reviewed in Sec. [31 where both a spherical basis, which is the direct extension of the Pauli 
matrices to three components, and a Cartesian basis with real components, are defined. The 
Maxwell equations are written in terms of the spherical spin matrices and combined into the 
Dirac equation form in Sec. HI In Sec. [5j transverse and longitudinal projection operators are 
defined and used to separate the Maxwell equations and solutions into the corresponding 
disjoint sectors. Lorentz invariance is addressed in Sec. [HI where transformations of the 
coordinates and derivatives, transformations of the Maxwell equation, and transformations of 
the solutions are explicitly written. In Sec. U\ plane-wave solutions, which are eigenfunctions 
of the momentum operator as well as the Hamiltonian, are given for both transverse and 
longitudinal states, and the set of solutions is shown to be complete. The explicit action of 
Lorentz transformations on the plane- wave solutions is described. Properties of normalizable 
wave packets formed from the plane-wave solutions are illustrated. The angular-momentum 
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operator and the corresponding eigenfunctions are given and shown to be complete in Sec. [HJ 
In Sec. El the Maxwell Green function is written as an integral over the plane- wave solutions 
in a form analogous to the Dirac Green function. As examples of applications of the Maxwell 
Green function, formulas for radiation from a point dipole source and from a Dirac current 
source are derived in Sec. [TU1 A summary of the main points of the paper is in Sec. [TT] and 
brief concluding remarks are made in Sec. [T2"l 

The relation of the present study to earlier work is indicated in the sections where the 
particular topics are discussed. 

2. Three-vector Maxwell equations and the Dirac equation 

The Maxwell equations in vacuum, in the International System of Units (SI), are 

VE = £ (1) 

VxE+— = 0, (3) 
VB = 0, (4) 

where E and B are the electric and magnetic fields, p and J are the charge and current 
densities, e and /i are the electric and magnetic constants, and c = (eo/io) 1 ^ 2 is the speed 
of light. The continuity equation 

l + V.J = (5) 

follows from Eqs. ([T]) and (j2J). 

The form of the Maxwell equations considered here is analogous to the Dirac equation 
for the electron. The Dirac wave function <f>(x) is a four- component column matrix that is 
a function of the four-vector x. For a free electron, the Dirac equation is 

(ih'fd ll -m e c)<l>(x)=0 1 (6) 

where h is the Planck constant divided by 2ir, m e is the mass of the electron, 7 M , \i = 0, 1, 2, 3, 
are the 4x4 Dirac gamma matrices, given by 

n i 1 Y „.i ( < j " 



1 = ' -I ; Y= -a< M = 1 ' 2 ' 3 ' ( 7 ) 



/ and are the 2x2 identity and zero matrices, a 1 , % = 1, 2, 3, are the Pauli spin matrices 

! / 1 \ 2 / -i \ 3 / 1 

° = { 1 o J' a = (, i o J' a = (,0 -1 
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and the derivatives <9 M are 

We take the metric tensor g^ u to be 

0°° = l ; fl-" = -l, i = l,2,3; <T = 0, /x^i/. (10) 
In terms of the spin matrices, the derivative term in Eq. can be written as 

= ( 7 & a T v d y (ii) 

The Pauli spin matrices act on two-component spin matrices in the electron wave func- 
tion. Oppenheimer has suggested that since the Maxwell equations involve three-vectors, 
three-component matrices should be considered for constructing a photon wave function j^] . 
Here we implement such an extension by replacing the Pauli spin matrices in the Dirac 
equation by the analogous 3x3 matrices described in the next section. 

3. Three-component spin matrices 

As is well known, three-vectors and operations among them are interchangeable with 
three-component matrices and matrix operations. In this section, formulas for these matrices 



relevant to subsequent work are given. Some of these formulas have been given in [14| . It is 
useful to define both Cartesian and spherical matrices to represent three-vectors. 
The Cartesian matrix representing a vector a may be written as 

a 1 

a 2 1 (12) 

a 3 

where a 1 , a 2 , a 3 are the rectangular components of the vector a, and a spherical represen- 
tation is denoted by 

a s = Ma c , (13) 

where M is a 3 x 3 unitary matrix specified in the following. The dot product of two vectors 
is 

a ■ b = a% = a%, (14) 

where f denotes the combined operations of matrix transposition and complex conjugation. 

Explicit Hermitian r matrices (r* = r), which are 3x3 versions of the Pauli matrices, 
are obtained by taking r 3 to be diagonal 

1 

r 3 = [ J (15) 
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and applying appropriate rotation matrices to obtain r 1 and t 2 : 

/ 1 

r 1 = 2 )( 1 )({0,f,0})r 3 2)( 1 )({0,-f,0}) = -= 10 1], (16) 

^ 2 \ 1 

. / -1 \ 

r 2 = 2)( 1 )({0,0,f})r 1 2)«({0,0,-f}) = ^= 1 -1 , (17) 

V2 \ 1 / 

where D^\{a, (3, 7}) is the j = 1 representation of the rotation group, parameterized by 
the Euler angles a, [3, 7 [lil ]. In particular, ©^(jO, | ,0}) represents the rotation about the 
2 axis by the angle ti/2 and ©^({O, 0, |}) represents the rotation about the 3 axis by the 
angle tt/2. The same rotations starting from a 3 , with the j = \ representation, reproduce 
a 1 and a 2 . The r matrices are related by 

[r\r j ] =\e l]k r\ (18) 

where e^k is the Levi-Civita symbol^] 

The cross product of two vectors a and 6 can be written in terms of the scalar product 
of the tau matrices with the vector a 

ra = tV (19) 

acting on the spherical matrix for the vector b as 

rab s = i (a x b) s , (20) 

provided the matrix M in Eq. ( TT3"j) is suitably chosen. To determine M, we take the 
Cartesian definition 

(0x6)* = e tjk a^b k (21) 

and write Eq. (1201 as 

raMb c = iM(axb) c . (22) 



1 

1 The tau matrices defined by Oppenheimer [5|] are N^t 1 AT, where N = [ i ] . [The minus sign 

v -1 

in r 2 in Eq. (10) of that paper apparently is a typographical error, as indicated by inspection of Eq. (11) 
The matrices defined by Majorana [l6| are M^t 1 M, where M is given in Eq. 
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Imposing the requirement that this equation be valid for any vectors a and b fixes M, up 
to a phase factor, to be 




M = — | V2 I , (23) 

which yields 

-^-ia 2 ) 

a s = I a 3 I . (24) 



^+ia 2 ) 



Consequences of Eq. (J2"0~j) are 



t ■ a, a s = 0, (25) 

r ■ a 6 S + t • 6 a s = 0, (26) 

alrbc s = ia-(bxc), (27) 
(r • a Cg) • (r • 6 d s ) = (a x c) • (6 x d) 

= abed— adcb. (28) 

Equation (|28l) can be written as 

4{t ■ o) ] t ■ bd s = c\(a-b-b s ai)d s (29) 

for any vectors c and d, which yields the relation 

(r • o) ] t b = a b- 6 s a s f , (30) 

where it is understood that the first term on the right includes the 3x3 identity matrix as 
a factor and the second term is also a 3 x 3 matrix. If a c has real components, then 

( T -a) ] = ra, (31) 
(r - a) 3 = a 2 t • a, (32) 

where a 2 = a ■ a is the ordinary real vector scalar product. 
Real Cartesian tau matrices r may be defined so that 

t ■ a b c = {axb) c . (33) 

This relation follows from Eqs. fl20l and fflBl with the definition 

t 1 = -iMVM, i = 1,2, 3. (34) 
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These matrices are antisymmetric t t = — f, where T denotes matrix transposition, in 
contrast to = r. For vectors a and b with real Cartesian components, we have 



t ■ a t ■ b 

\3 



b c a c 



r ■ a 



[t ■ a 



a ■ b, 

a 2 (f ■ a), 



(35) 
(36) 
(37) 



The matrix 



r ■ a 






-a 3 


a 2 


a 3 





—a 


-a 2 


a 1 






(38) 



has the form of the lower right portion of the electromagnetic field-strength tensor F^ v ', as 



given in [17j for example. 



4. Matrix Maxwell equation 

In terms of the notation of the previous section, the matrix forms of the Maxwell equa- 
tions in d2J) and ([2D, for the source-free case (J = 0), are 



iT-VB s + 



c dct 

1 T ■ V^s - C- 



dct 



0. 



0. 



These equations may be written as two uncoupled equations 

d 



dct 

A 
dct 



t- V 1 [E s + icB s 
t- VI (£ s -iaB s 



0. 



(39) 
(40) 

(41) 
(42) 



where J is the 3x3 identity matrix. In the Cartesian basis, Eqs. (jUJ) and (J52J) are 

I ^ct +if ' V ) (^ + icfi c) = 0, 
. v] (£ n -iaB n ) = 0. 



(43) 
(44) 



In these expressions, it is evident that for real electric and magnetic fields, Eqs. (143]) and 
(jSJ) are complex conjugates of each other and reduce to a single complex equation. It was 
recognized in lectures by Riemann in the nineteenth century that this complex combination 
of E and B is a solution of a single equation 18]. This fact was also discussed in 19, 2(| and 

is included in many works up to the present. Equations (143]) and (144]) may be interpreted as 
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Maxwell equations for right- and left-circularly polarized radiation, analogous to the Weyl 



equations for right- and left-handed neutrino fields [8|, [14 . 

However, in this paper, we consider the more restrictive case of complex electric and 
magnetic fields that are simultaneously solutions of both Eqs ( ]4"5|) and (T44|) . or equivalently 
both Eqs. (I39p and (1401) . for any polarization of radiation. The question of whether such 
solutions can be found is answered by their explicit construction in subsequent sections of 
the paper. To formulate this approach, we follow the Dirac equation and write 

'& -*VM = o, (45) 

which is a restatement of Eqs. (1391) and (1401) in the form of the Dirac equation for an electron 
wave function. It is a matrix equation with six components that may be viewed as a single 
equation equivalent to Eqs. (1391) and (T40]) for any polarization of the fields. Any complex 
solution of Eq. (T4"5])_ is a solution of both Eqs. (T4T1) and (T4"2l) . Similar wave functions have 
been discussed in [2l|, |22|, |23j . It should be noted that this formulation is different from the 
six-component form considered by Oppenheimer in which the upper-three components and 
lower-three components represent opposite helicity states 
If we define 6x6 gamma matrices by 

V=(J J): 7^(_° T ), -1,2,3, (46) 
where is the 3x3 zero matrix, and write 




then Eq. (145]) takes the covariant Dirac equation form 

-fd^ix) = 0, (48) 

which provides a concise expression for two of the Maxwell equations. We can also write 
this as 

V{xfd^ = 0, (49) 

where W(x) = &^(x) 7 and d M denotes differentiation of the function to the left. Although 
these equations are simply algebraic rearrangements of the two Maxwell equations, the 
resemblance to the Dirac equation and wave function is suggestive of a form that photon 
wave functions might take. 

It is of interest to note that for solutions of the Dirac equation for the hydrogen atom, 
the lower two components are small and approach zero in the nonrelativistic limit, i.e., as 
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the velocity of the bound electron approaches zero. Similarly, for local electromagnetic fields 
generated by moving charges, the magnetic field, given by the lower three components of 1^, 
also approaches zero in the limit as the velocity of the charges approaches zero. 
To take source currents into account, Eq. (|2J) is written as 

1 BW 

ir-VB 8 + --^ = -noJ B , (50) 
and a source term S is defined to be 

*(») = ( -" 0C J ' W ) . («) 

where is a 3 x 1 matrix of zeros. This yields the expressions 

-fd^{x) = E{x) (52) 

and 

Vixfdrf = S(x), (53) 



either of which is referred to as the Maxwell equation here. The source term in Eq. (jo2|) 
or (|53p represents a fundamental difference between the Dirac equation and the Maxwell 
equation, as mentioned in Sec. Q] (5j]. 

In this framework, an energy-momentum density operator is 



p" = |V, (54) 



which gives 



.0 ,r, _ i,i rr-12 1 



Vcp u V = - [e \E\ 2 + — \B\ 2 ) =u, (55) 
2 V 



i£o 

2 
1 

C 2 /i 



ReExF = g. (56) 



Eqs. fl!3]) and ([53]) imply that 

B^(x)Y &(x) = E(x)&(x) +~&(x) S(x), (57) 
which is a complex form of the Poynting theorem [see Eq. ([HP ] 

— + V • S = -Re £ • J, (58) 
at 

where 

5 = c 2 g, (59) 



which gives the conventional result if the fields and current are real [17 
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5. Transverse and longitudinal fields 



To make a Helmholtz decomposition of electromagnetic fields expressed in matrix form 
into transverse and longitudinal components, we define 3x3 matrix transverse and longi- 
tudinal Hermitian projection operators iTj(a) and 77^ (a) to be 



nj(a) = ' T (60) 
JlJ-(a) = ^ (61) 



Based on identities in Sec. El these operators have the following properties: 

r 

i 2 



[nJ(a)Y = iTj(a), (62) 

[iT^(a)] 2 = iT s L (a), (63) 

JT s T (a) + i7 s L (a) = I, (64) 

JTj(a) iI s L (a) = 0, (65) 

nj(a)a s = 0, (66) 

n^(a)a s = a s . (67) 

Acting on the matrix of an arbitrary vector b, the operators project the components per- 
pendicular to and parallel to the argument a 

nj(a)b s = b s -—a s , (68) 

CI ' (X 

n^(a)b s = ^a s . (69) 

In addition to these algebraic relations, usefulness of the projection operators arises 
from an extension to include differential and integral operations acting on coordinate-space 
functions. Formally, we write 



i7 s T (V) = (70) 
iT^(V) = (71) 



which takes into account that fact that V has real Cartesian components (in the sense that 
they give real values when acting on a real function). The inverse Laplacian is defined by 
the relation 



which yields 



4/0*0 = -j- fdx'—±—f(x% (72) 
V z Air I \x — x \ 



V 2 ^/(x) = -^/d-'V 2 ^- f(x') = f(x), (73) 
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based on 



\x 



X' 



-Atx 6(x 



X 



(74) 



where 



S(x-x') = Six 1 -x n )S(x 2 -x ,2 )S(x 3 -x' 3 ). 



(75) 



Equation (I73p indicates that the Laplacian operator follows analogs of the rules of algebra 
in this context. For example, we have 



• (r-V) : 
V 2 

T ■ V) 2 



V 2 (r ■ V) 
V 2 V 2 



n (v) 



(76) 



where either simply canceling the V 2 factors in the numerator and denominator or applying 
the definition in Eq. (1721) for the operators acting on a suitable function gives the same 
result. Transverse and longitudinal components of the electric and magnetic fields and the 
current density are identified by writing 



where 



if + F S L , 



(77) 



(78) 
(79) 



and F s may be any of E s , B s , or «7 S . 

The separation of the Maxwell equations into two independent sets of equations which 
involve either transverse components or longitudinal components takes the following form. 
In terms of the spherical matrices, Eq. ([!]) is 



eo' 



10) 



where the transverse component of the electric field is absent, because VjiTj(V) = 0. 
Equations ([!]) and flHUj) are equivalent, in the sense that each can be derived from the 
other; Eq. ([T]) follows from Eq. (18"U|) if the vanishing transverse component is added to 
the latter equation. In the separated form, it is evident that the equation neither contains 
information about or places any constraint on the transverse component Ej. The transverse 
and longitudinal projection operators acting on Eq. (150|) . the matrix form of Eq. ([2]), yield 



ir-VBj + 



ldE2 



c dct 
IdE]^ 



T 

s i 



(81) 
(82) 



respectively, which take into account the commutation relation [iTj(V), r • V] =0 and the 
fact that JJg (V) t • V = 0. Together, these equations are equivalent to Eq. (1501 which can 
be restored by writing the sum of Eq. (1HT1) and Eq. (182]) and adding the term that vanishes. 
Evidently, this pair of equations is independent of B^. Similarly, Eq. ([3]), or equivalently 
Eq. (l4"Uj) . can be written as the pair 

ir-V£ s T -c^ = 0, (83) 



Oct 



0, (84) 



Oct 

which are independent of E^. Equation (J4]) takes the form 

V|B S L = 0, (85) 

independent of Bj. The transverse and longitudinal equations comprise two independent 
sets. 

Six-dimensional transverse and longitudinal projection operators are defined by 

iTJ(V) 



tf T (V) = „ ' , (86) 



o iijrv 



" s 



il^(V) o 



n L W = ^ , (87) 



i7 = L (V 



' s 




where is the 3x3 matrix of zeros, i7 T (V) + 77 L (V) = J, and J is the 6 x 6 identity 
matrix. The transverse equations are summarized by writing 



Yd^ T (x) = ~ T (x), (88) 

where 

V T (x) = i7 T (V)!f(x), (89) 

~ T (x) = 77 T (V)~(x). (90) 

Equation (I8"8"|) also follows directly from Eq. (152"|) and the fact that [77 T (V), 7 M 9 M ] = 0. 
The longitudinal equations are Eqs. (f8"0l . (18"2"|) . and (IHol) . together with the continuity 
equation, Eq. ([5]), which can be expressed as 

% + ViJ s L = 0. (91) 

Since the continuity equation follows from Eqs. flHUl) and (I8"2l . it is not necessary to include 
it in an independent set of equations; it is listed here only to show that it provides no 
restriction on «7 S T . Equations (I8"l"|) and (I8"5"|) are eliminated from consideration by taking 



(92) 



Constant fields are eliminated by the requirement that static fields vanish at infinite distances 
for finite source distributions. However, a constant magnetic field may be approximated by 
the field at the center of a current loop with a radius that is large compared to the extent 
of the region of interest. Such a steady-state current density is transverse, as shown by 
Eq. (1821) . and so the magnetic field, given by Eq. (1ST]) , is also transverse, which is consistent 
with Eq. (|9"2l) . A complete set of equations, equivalent to the set of Maxwell equations, is 
provided by Eqs. fl52|) . flHUj) . and fl92l) . and the transverse fields are completely described by 
Eq. (|88|). 



6. Lorentz transformations 

Lorentz transformations of the matrix Maxwell equation are examined here in order 
to confirm that this form of the Maxwell equations is Lorentz invariant. We adopt the 
convention that transformations apply to the physical system rather than to the observer's 
coordinates. 

To represent four-vector coordinates, the Cartesian matrices are extended to include a 
time component x° = ct, so coordinate vectors take the form 

/ x° \ 

' ct 



X 



x 1 



x 2 



(93) 



V* 3 / 

We employ the Cartesian basis for coordinate and momentum vectors and the spherical 
basis for fields and currents, with a few exceptions that will be apparent. The notation x 
represents either the four-coordinate argument of a function or a column matrix, depend- 
ing on the context. It is sufficient for our purpose to consider only homogeneous Lorentz 
transformations and to consider rotations and velocity transformations separately. These 
transformations acting on four-vectors leave the scalar product 

x-x = x T g x = (ct) 2 - x 2 (94) 

invariant, where g is the metric tensor given by 

A remark on notation is that a boldface means either a 3 x 3, a 1 x 3, or a 3 x 1 
rectangular array of zeros, as appropriate. We take the liberty of using an ordinary zero on 
the right-hand side of equations to mean whatever sort of zero matches the left-hand side. 

6.1. Rotation of coordinates 

Rotations are parameterized by a vector u = 9u, where u is a unit vector in the direction 
of the axis of the rotation and 9 is the angle of rotation. An infinitesimal rotation 59 u 
changes the point at position x to the point at position x', where 



x' = x + 59 u x x + . . . , 
14 



(96) 



or 

x' c = (1 + 59? -u)x c + ... . (97) 

For a finite rotation, the operation is exponentiated to give 

x' c = e™x c = R c (u) x c . (98) 

Expansion of the exponential function in powers of 9, taking into account the fact that 
(r ■ u) 3 = —t ■ u, yields 

R c (u) = I + (T-u) 2 (l-cos9)+T-usin9 

= u c uj — (t ■ u) 2 cos 6 + f ■ u sin 9. (99) 

Evidently, i?~ 1 (w) = R c (—u) = Rj(u). It is confirmed that this operator has the appro- 
priate action on a vector by calculating 

x' = uu ■ x — u x (u x x) cos9 + u x xsm.9. (100) 

We use the notation 

x' = R(u)x (101) 

to represent the transformation in Eq. fllOOp . 

Rotations of a four-vector only change the spatial coordinates and are written as 

x' = R{u)x= I Ct I , (102) 



RJu) x ( . 



where 



Riu) ' o *„(„))■ (103) 

The scalar product x ■ x is invariant under rotations, since x 2 is invariant. 

The spatial coordinate rotation operator in the spherical basis, which follows from 

R s (u) = MR c (u)M\ (104) 

is 

Rs(u) = e~ [T U = u s ul + (t • u) 2 cos 9 - i r • u sin 9, (105) 

and -R ( 7 1 (u) = R s (—u) = Rl(u). Starting from the geometrical constraint that the rotated 
cross product of two vectors is the cross product of the rotated vectors, written as 

R s (u)(axb) s = {a'xb%, (106) 
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we have 

R s (u)t ab s = t ■ a' b' s = t ■ a' R s (u) b s . (107) 
Since this relation holds for any vector b, it yields 

Rs(u) r a = t ■ a'R s (u) (108) 

and 

R s (u)t ■ aR;\u) = t o!. (109) 

A direct calculation provides the same result. 

The relation between the rotated and the unrotated gradient operators is given by 

„,i d dx j d dx j , . « 

V = r= r = 110) 

dx' 1 dx' 1 dxj dx n v ; 



R-Uu) = Rcijiu), (111) 



= i* c (u)V c , (112) 



where, from Eq. (1981 . we have 

dxi 
~dx Ti 

so that 



and from Eq. (JTOlj) . 

= J?s(ti)V 8 . (113) 

Since the spherical gradient operator transforms as a spherical vector, we also have 

Rs(u)t- VR-\u) = t-V (114) 

from Eq. fllQ9p . Equations (11131) and (jl 14j) imply that transverse and longitudinal projection 
operators transform according to 

JTj(V') = R s {u)n^{V)R-\u), (115) 
iT s L (V) = R s (u)Il]:(V)R- s \u). (116) 

The action of the inverse Laplacian in terms of the rotated coordinates is the same as it is 
for unrotated coordinates, which follows either because V' 2 = V 2 from Eq. (11131) or by the 
definition in Eq. (1721) . taking into account the fact that the Jacobian for a rotation is unity. 
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6.2. Velocity transformation of coordinates 

Velocity transformations are parameterized by a velocity vector v = c t&nh(v. If a 
space-time point is given an infinitesimal velocity boost of S(cv, its spatial coordinate will 
change to 



x 



x + S( ctv + . . . , 



(117) 



and its time coordinate must transform in such a way that the scalar product is invariant. 
In particular, we require x' ■ x' = x ■ x, which yields 



ct' = ct + S( v • x + . . . . 
The complete infinitesimal transformation is 

' ct 



ct' 
x' 



1 Vc 



+ 



(118) 



(119) 



This may be written in terms of a 4 x 4 matrix valued function of the velocity direction: 

K(«) = ("/'), (120) 

for which 

* 2 (*' - ( o ) (121) 

and K 3 (v) = K(v). For a finite velocity, the transformation is exponentiated to give 

x ' = x = V(v) x. (122) 

Expansion in powers of ( yields 

V(y) = I + K 2 (v)(cosh( - 1) + K(v)sinh( 

cosh ( vj sinh ( 

v c sinh ( I + v c vj (cosh ( — 1) 



(123) 



The relations V T (v) = V(v) and gV(v) = V 1 (v)g confirm the invariance of the scalar 
product: 



x ■ x = x T V T (v) g V(v) x = x ■ x. 



The transformation yields 



ct 1 = ct cosh ( + v ■ x sinh (, 

x' = x + vv ■ a:(coshC — 1) + ctv sinhC- 
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(124) 



(125) 
(126) 



A point with x = has the boosted velocity 

x' 

— = ctanh(v = v. (127) 

t 

The spherical counterpart of the operator V(v), in the velocity transformation 

( i ) = *««> ( 1 • (128) 

is 

cosh £ -Oj sinh £ 

6 S sinh £ J + t) s v| (cosh ( — 1) 

For the four-gradient operator, we have 



(129) 



_d_ _ dx^_d_ = dx^ 

and from Eq. (11 22p . which can be written as 

x = V~\v)x' = V(-v)x' (131) 

or 

x v = V VjJb (—v) x" 1 , (132) 

we also have 

= V^(-v) = V^-v), (133) 

which yields 

% = V, u {-v)d u . (134) 
If a Cartesian gradient operator is defined as 

d c = ( ^ j , (135) 

then Eq. (11341) gives 

<7# c = V(-t0<7d c (136) 

or 

# c = V(v)d c , (137) 

since (7 V(—v) g = V(v). 
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6. 3. Parity and time reversal of coordinates 

Lorentz transformations that leave the scalar product in Eq. (1941) invariant include the 
parity transformation P = g, time reversal T = —g, and total inversion PT = —I operations. 
These transformations have the following defining effects on the coordinate vectors: 

Pr = : V (138) 



X 



T X = ( ;f ) , (139) 
PTx = -x. (140) 

It is sufficient for the present purpose to consider only P and T. The coordinate derivatives 
transform as 

Pd c = f ^* Y ( 141 ) 
Td c = ^ ~~$ct \ . (142) 

Comparison of Eqs. fl95l) and (I103P shows that parity transformations commute with 
rotations. On the other hand, for velocity transformations, the relation 

PV(v) = V(-v)P (143) 

applies as it should, because the space reflection of a point moving with a velocity v is a 
point at the reflected position moving with a velocity — v. Similar conclusions follow for 
time-reversal transformations. 



6.4. Rotation of ]P(x) 

The result of a rotation, parameterized by the vector u, applied to the field \P(x) in 
Eq. f)52p is the field W{x) given by 

&{x) = Tl{u)^{R- l {u)x), (144) 

where lZ(u) is a 6 x 6 matrix that gives the local transformation of the field \P(x) at any point 
x. The inverse transformation of the argument on the right-hand-side takes into account the 
fact that the transformed field at the point x originated from the field at the point that is 
mapped into x by the transformation. Lorentz invariance is confirmed by showing that the 
transformed field satisfies the same equation as the original field. We expect the current to 
transform in the same way as \P and write 

E'(x) = TZ(u)S(R-\u)x). (145) 



19 



The objective is to show that 

-fdp#\x) = E'(x), (146) 

for a suitable transformation lZ(u). In terms of the original field and source, Eq. (11461) is 
given by 

7^^(w)^(i?" 1 (n)x) = K{u) ~(R-\u)x) (147) 

or 

Y$Jl(u)$(x) = K(u)E(x), (148) 
where the variable x has been replaced by x' = R(u) x. Thus Eq. (11461) will follow if 

n-\u)^d'Jl(u) = Yd,. (149) 
We expect lZ(u) to be of the form 

*<«> = ( T RJu) ) ■ ^ 



which yields 



iJ-r RlH^T -V'RJu 



so Eq. (1 149ft follows from 

R^i^r -VR s (u) = r-V, (152) 

which, in turn, follows from Eq. (11141) . We conclude that, as expected, the solution and 
source terms, transformed according to Eqs. (11441) and (I145j) . where TZ(u) is given in 
Eq. (I150p . satisfy the same equation as the original solution and source terms. The six- 
dimensional rotation operator lZ(u) may be written as 

TZ(u) = e~ lSu , (153) 

where 

';:). (164) 

Equations (11441) . (11451) . and (I150p correspond to the separate equations 

E' s (x) = RbME^R-^x), (155) 

B' B {x) = R B (u)B e (R- 1 (u)x), (156) 

J' s (x) = R s (u)J s (R- l (u)x). (157) 
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It can be confirmed that Eqs. ([T]) and (jlj) in spherical form are invariant under rotations. In 
particular, Eq. (pQ) for the rotated electric field and charge density p'(x) = p(R~ 1 (u)x) is 

V|^(x) = 

or 

VtRME^R-^u) x) 
The substitution x — > R(u) x gives 

V? Rs(u)E 8 (x) 

and Eq. (jTT5j) yields 

VtE s (x) = 

Hence, the transformed field and charge density satisfy Eq. ([I]) if the original field and charge 
density do. Similarly, 

VtB' s (x) = VtB s (x) = 0. (162) 

Thus, all of the Maxwell equations in matrix form are invariant under rotations. 

The separation into transverse and longitudinal components of the electric and magnetic 
fields is also invariant under rotations. This can be seen by considering the expression 
iT s (V) F s (x), where iT s (V) is either i7j(V) or Il\ (V) and F s ( x) is any of E s (x), B s (x), 
or J s (x). We have 

27 s (V)F s '(x) = nWRtMF^R-^x) (163) 

or 

n s (V')F^x') = n s (V')R s (u)F s (x) 

= R s (u)n s (V)F s (x), (164) 

where the last line follows from either Eq. (I115P or Eq. (11161) . This means that if the 
original field is transverse or longitudinal, then the rotated field has the same character. 
These results extend directly to the six-dimensional projection operators i7(V), solution 
lP(x), and source 3{x). 



P'{x) 



(158) 



p(R- l (u)x) 
eo 



(159) 



eo 



(160) 



p(x) 
eo 



(161) 
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6.5. Velocity transformation of \P(x) 

The result of the velocity transformation, by a velocity v, applied to the field \P(x) in 
Eq. (|52|) is the function $'(x) given by 

$'{x) = V(v)^V-\v)x), (165) 

where V(u) is a 6 x 6 matrix that gives the local transformation of the field \P(x) at any 
point. The inverse transformation of the argument on the right-hand-side plays the same 
role as for rotations. Our objective is to establish the covariance of Eq. (1521) by showing that 
if ty{x) is a solution of that equation with a source S(x), then 

rd^'(x) = E'(x), (166) 

where E'(x) is a suitably transformed source term. Equation (11661) can be written as 

7^V(v)^(V" 1 (v)x) = E'{x) (167) 

or 

7^V(v) <%) = E'(V(v)x), (168) 

where the variable x has been replaced by x' = V(v) x. The 6x6 matrix V(v) is based on 
the conventional local velocity transformation of the electric and magnetic fields as discussed 
in [A] Here we write it as 

V(v) = e iK *, (169) 

where 

K=(l nV (170) 



V T , 

Expansion of the exponential function in Eq. (I169p in powers of ( yields 

V(v) = J+(/C-v) 2 (coshC-l) + (/C--0) sinhC 

(r ■ v) 2 (cosh ( — 1) t ■ v sinh ( 

t ■ v sinh ( I + (t ■ v) (cosh £ — 1) 

where (K ■ v) 3 = K. ■ v is taken into account, so that 



, E' \ ( ^ s + (r-i)) 2 J E; s (coshC-l) + ir-'0cS s sinhC 
V(v)$ = = , • (172) 

1 icB' s \ i[cB s + (T-v) 2 c J B s (coshC-l)-iT-v J E; s sinhC] 



In Eq. (11681) . we have 

J Set' T V 
22 ac 



7^; = ( dct^ ' y Q ] ; (173) 



where, from Eqs. (I135p and fl!37j> . 



d d 

cosh£- sinh^'D-V, (174) 



dct' dct 

d 

V = V + (coshC - 1)* v ■ V - sinh (6 t— . (175) 

act 

Multiplication of Eq. ffTTTj) by Eq. (jT73j) yields the identity 

7 °» V{ } V -T -|- "0 s t)| (cosh — 1) J 



— sinh ( v s ^i 

sinhC^Vl 



and hence 



T»ffV(v)9(x) = ( / + ^ s t (coshC-l) . 

' " I / + v 8 vl (cosh C - 1) I K ' 



+ 



— sinh ( v s V • -E(x) 
i sinh £ v s c V • -B(x) 



-/i c [J s (x) + v s v ■ J(x) (cosh ( — 1) + sinh ( v s cp(x)] 





-p Q cJ' B {x) 




where J' s {x) is the velocity transformed three- vector source current, and 



(176) 



(177) 



S'(x) = (-^™*)\ (1?8) 

The result in Eq. (11771) takes into account the additional two Maxwell equations in (pQ) and 
(jlj), besides Eqs. <^fy and ^ used to construct Eq. (1521) . It also requires the conventional 
result that the current 



cp(x) 
J B (x) 



(179) 



transforms as a four- vector under the velocity transformation given by Eq. (11291) . Equa- 
tion (I177p establishes the validity of Eq. (I166p . provided the transformed three- vector source 
current is the three- vector component of the transformed four- vector current. 

The covariance of Eq. fl52l) , even though the charge density does not appear in the source 
term, is linked to the fact that the current density satisfies the continuity equation. Since 
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the continuity equation follows from the Maxwell equations, it cannot be expected that 
consistent solutions may be found for an arbitrary four-vector current density. However, for 
valid sources, information about the charge density may be obtained from the three-vector 
current density and the continuity equation. For example, if the electric field is specified 
at a particular time, then the charge density at that time is known from Eq. ([T]) and may 
be determined at any other time from knowledge of the three-vector current density by 
use of the continuity equation. Thus the time evolution of the electromagnetic fields can 
be described relativistically with no reference to the charge density. The Maxwell Green 
function, which provides the solutions of Eq. (152]) . is discussed in Sec. [9j 
Since V^7°V = 7 and V^r/V = j°rj, where 



the invariance of the quantities 



'/=( °j )" 11805 



W = \E\ 2 -c 2 \B\ 2 , (181) 
~Wr)& = 2icReE ■ B (182) 



is evident. 



6.6. Parity and time-reversal transformations of \P(x) 

We expect the fields ty{x) to transform under a parity change according to 

$\ x ) = V^P^x), (183) 

where V is a 6 x 6 matrix, and we assume that the current transforms in the same way, so 
that 

S'{x) = VS{P~ 1 x). (184) 

We can obtain 

-fd^\x) = E'(x), (185) 

by finding a suitable matrix V . In terms of the original field and source, Eq. (1185P is given 
by 

Yd^WiP^x) = VSiP^x) (186) 

or 

Y&^V&ix) = VS(x), (187) 

where d'^ is the parity transformed derivative given by Eq. (11411) . Equation (I185p will follow 
if 
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(188) 



which corresponds to 



V' x iV = -f, 1 = 1,2,3. (190) 
Solutions of these equations are provided by 

V = ±(J _° J ). (191) 

The minus sign corresponds to the conventional choice of how classical electric and 
magnetic fields transform under a parity change, that is, the current and electric fields 
change sign and the magnetic fields do not (l7| . 

To examine time-reversal invariance, we first consider field which is real in 

the Cartesian basis. In this case, the conventional use of an anti-unitary operator for time 
reversal is unnecessary, and the same can be expected to be true for fields expressed in the 
spherical basis. We write 

V\x) = T^T^x) (192) 

and 

E'{x) = —TE(T~ l x), (193) 

where T is a suitable 6x6 matrix. The minus sign for the current provides the result 
that the electric field does not change sign under time reversal and the current does. The 
objective is to find a matrix T such that 

Yd^\x) = S'(x) (194) 

or 

^df/T W{T~ x x) = -T~(T-V), (195) 

and so 

Y^T^ix) = -TS(x), (196) 
where d'^ is given by Eq. (11421) . Such a matrix satisfies 

T-y^T = -Yd, (197) 

or 

T-yT = 7°, (198) 

T-yr = - 7 \ i = 1,2,3, (199) 
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with a solution provided by 



T = (o -/)• < 200 > 

The matrices V and T commute with the rotation matrix, as they should, and we have 
the result that the matrix form of the Maxwell equations is invariant under parity and 
time-reversal transformations. 

For quantum-mechanical time reversal, the time-reversal operator is anti-unitary and 
includes complex conjugation, or Hermitian conjugation in the case of a matrix solution, 
which has the effect of interchanging initial and final states. For an example where such an 



interchange corresponds to observable consequences in QED, see |24l. |25| . We thus write 



W'{x) = %^{T~ x x) = y ] (T- l x)U~ l , (201) 
= -l-iT^x) = -~\T- l x)U~\ (202) 



E'(x) 



where T = CU is the product of the Hermitian conjugation operator C, which has the action 
C &(x) = &^(x) and a unitary matrix U. The objective is to find a U such that 

V'(x)7°*d^ = ~'(x)7° (203) 

if \P(x) is a solution of Eq. (j55|) . Equation f)203p can be written as 

!F t (a?)W- 1 7°l?Vy' 4 = -E\x)U~ 1 1 , (204) 
where d'^ is given by Eq. (11421) . Equation (12041) follows from Eq. (J53"j) provided 

^W^y^U^ = -^ t7 M , (205) 

which has as a solution 

U = T (206) 

and 

% = CT. (207) 

7. Plane-wave eigenfunctions 

Following the analogy with the Dirac equation, the Hamiltonian for the Maxwell equation 



is 



H = cctp = -ihcct-V, (208) 

where a 1 = 7° 7*, and the wave functions for the photon may be identified as the complete 

set of eigenfunctions of 7i. The solutions considered here are coordinate-space plane waves 
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characterized by a wave vector k and a polarization vector e A ; both positive- and negative- 
energy solutions, as well as zero-energy solutions are included to form a complete set. These 
solutions are also eigenfunctions of the momentum operator 

V = lp = -ihlV, (209) 

which commutes with H.. 

The plane-wave solutions are not normalizable, because their modulus squared is inde- 
pendent of x and the integral over all space does not exist. As a result, the solutions include 
an arbitrary multiplicative factor, that could be a function of k. Here a factor is chosen to 
provide the simple result in Eq. (12251) . 

7.1. Transverse plane-wave photons 

We first consider transverse photons, i.,e., photons for which the electric and magnetic 
fields are perpendicular to the wave vector. The polarization vector is a unit vector pro- 
portional to the electric or magnetic fields, represented by a three component, possibly 
complex, vector in the spherical basis. As such, the polarization vector does not transform 
as the spatial component of a four-vector under velocity transformations. 

Two polarization vectors, both in the plane perpendicular to k, are denoted by 

e A (fc); A = 1,2 

They have the orthonormality properties 

e\ 2 (k)e Xl (k) = 5 X2M , 
ftfo(fc) = 0, 

and the completeness property 

2 

^e A (fe)4(fe) = I-k s kl = (r-k) 2 = nJ(k) 

A=l 

From Eq. (12 12p . we also have 

(r-fc) 2 e A (fc) = e A (fc). (214) 

The polarization vectors can represent linear polarization, circular polarization, or any com- 
bination by a suitable choice of e x (k). For example, for k in the e 3 direction, linear polar- 
ization vectors in the e 1 and e 2 directions are 




(215) 
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(210) 



(211) 
(212) 



(213) 



according to Eq. (|2^|) . Similarly, circular polarization vectors are (see Sec. 18. 2j) 

ei(e 3 ) = | ] ; e 2 (e 3 ) =| I. (216) 





These polarization vectors can be transformed to the vectors corresponding to any direction 
of k with the rotation operator in Eq. (11051) . (See also Sec. 17.51 ) 

Positive (+) and negative (— ) energy transverse photon wave functions are given by 

fc,AV ; V^) 3 V r-k e x {k) J v ; 

Although constructed geometrically to be transverse by the choice of polarization vectors, 
these wave functions are also transverse in the sense defined in Sec. [5j We have (see [B] for 
more detail) 

n T (v)^ x (x) = n T (k)^l(x)=^ x (x), (218) 
n L (v)^ x (x) = n^k)^ k ± l(x) = o, (219) 

where 77 T and IJ L are defined in Eqs. (IBS"]) and (IH71) . The wave functions in Eq. (I217p are 
eigenf unctions of the Hamiltonian in Eq. (12081) with eigenvalues ±hc \ k\. In particular, 



Wi3(a0 = ±hc a -k^(x), (220) 



and 



r-k\f e x (k) \ _ |fc| / e x (k) . (221) 



so that 



Also, 



r ■ k ) \ T-k e x {k) ) 1 1 V r • k e x (k) 

H^l(x) = ±Hc\k\^(x). (222) 



T^ x \x) = ±hk^(x). (223) 
The wave functions have the expected property 

^toO*) = 0. ( 224 ) 

since the electric and magnetic field strengths are equal for a transverse photon. Normal- 
ization and orthogonality relations are 

jdx (aOVffk (*) = W(*2-*i), (225) 

/ dx^l(x)^ Xl (x) = 0. (226) 
7 28 



The latter relation follows from a cancellation of terms between the upper-three and lower- 
three components of the wave function: 



'.(=F) 



J + r ■ k 2 r • fci 



=Fi(fc 2 +fci)-a; 



J + r • fc 2 r • ki 



e Al (fci)5(fc 2 + fci) = 0. (227) 



The transverse wave functions constitute a complete set of such functions. The com- 
pleteness is established by writing 



J2[dk C^CVi) 

A=l J 



A=l 



E/ dfc 



e A (fc) gT(fc) 



i x (k) e{(k)r ■ k \ ± ik .( X2 - Xl) 



2(2vr) 3 V r ■ k e x (k) e\(k) r ■ k e x (k) e\(k)r ■ k 



(t ■ kf 

t ■ k 



t ■ k 

(T-fc) 5 



_,±ifc-(cc2 — £Ci) 



(228) 



To evaluate the integrals for the sum over positive and negative energy solutions, we use k 
to represent either a plus sign or a minus sign and write 



and 



[ dk(r- k) 2 e Kik < X2 - x ^ = 2(2tt) 3 iT s T ( V 2 ) S(x 2 - Xl ) 
E / dkr-ke Kik < X2 - Xl) = 0, 



(229) 



(230) 



which yields the transverse completeness relation 



E E / dk <i(*2)<f (*i) = ^ T (v 2 ) 5( X2 - Xl ) 



(231) 



K->± A = l 



7.£. Longitudinal plane-wave photons 

The transverse wave functions alone do not provide a complete description of electro- 
magnetic fields. For example, the field of a point charge q at rest at the origin, given by 



E s {x) 

is purely longitudinal, because 



q x s 



47re n \x\ 



q -V -1- 

S I I ) 



47re kc 



IT B L (V)^ B (x) = 

ilJ(V) E s (x) = 0. 
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(232) 



(233) 
(234) 



The longitudinal photons are represented by a third polarization state, labeled A = 0, with 
the polarization vector taken to be 

6o(fe) = fe s . (235) 
If k is in the e 3 direction, the longitudinal polarization vector is (up to a phase factor) 



eo(e 3 ) = 1 . (236) 




Longitudinal wave functions are 

*" {x) = 7?b ( e f' ) e ' k " (237) 



or 



■4>tt{x) = " f )e~ ikx . (238) 

This polarization state has the property that 

n\v)^{x) = n\k)^{x) = ^{x\ (239) 

n T (V)^(x) = n T (k)^(x) = 0. (240) 
The wave function has an energy eigenvalue of zero, 

H^(x) = ±tica.kil>W{x) = 0, (241) 
because r • ke (k) = 0. However, 

V^(x) = ±hk^(x). (242) 

Normalization and orthogonality of the A = wave functions, as well as the transverse 
wave functions, are given by Eqs. (12251) and (j226| . where Ai and A2 may take on any of the 
values 0, 1, or 2. 

Completeness relations are given by 

dfc = Jdk JL ( l ) e-(— ) ( 2 43) 



and 



/<* *kWfc,»(*0 = /<* ^ ( ° . o( ^°.t (S) ) .-*<—>, (244) 



where 



which yields 



J dk e (k) e j {k)e ±ik < X2 - x ^ = (2tt) 3 i7^(V 2 ) 5(x 2 - Xl ), (245) 
J dk 4%x 2 )^( Xl ) = i7 L (V 2 ) 5(x 2 - x x ). (246) 



For the example of a point charge at the origin, we have 



*p(*) = 7— IZiifl n S )> (247) 



47re |cc| 3 \ 
which may be written as (see O for some detail) 

& p (x) = n L (V)V p (x) = /da:iI7 L (V)<?(x-Si)!F p (a;i) 
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J^w\^ x) - (248) 



7.5. Fit// orthogonality and completeness of the plane wave solutions 
The full orthogonality relations are 

Jdx^ 2 (x)^l(x) = <W<W(fc 2 -fci), (249) 

where the factor 5 K2Kl is 1 if n 2 and K\ represent the same sign and is for opposite signs, 
and A 2 and Ai may be any of 0,1,2. The combined result of the transverse and longitudinal 
completeness relations, Eqs. (I23ip and (12461) . is 

EE/ dfe W = I6 (^ - *0> (250) 

A=0 

where 77 T (V) + i7 L (V) = J. 

7.^. Time dependence of the wave functions 

The time dependence of the photon wave functions is given b}{§ 

V£i(*) = V>£(*) e^', (251) 



2 The notation /(a?) = /(x)| f _ is employed throughout this paper. 
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where u is determined by the equation 



rd^ k K l(x) = /(j^ + a-vj^HO. (252) 
For transverse photons 

u = c\k\, (253) 

and for longitudinal photons 

uj = 0. (254) 

The complete exponential factor is thus 

e~ Kik - x , (255) 

where 

k - ( t ) « ( I) • < 256 > 

depending on whether the photon is transverse or longitudinal. This corresponds to the 
eigenvalue equation 

= khu>4>%{x) (257) 

and a time dependence given by 

= e-^ViSC*). (258) 

It is also of interest to consider the effect of a hypothetical photon mass m 7 on the 
longitudinal photon wave function. Such a modification, with an infinitesimal mass, resolves 
an ambiguity in the construction of the Green function as discussed in Sec. M Following the 
form of the Dirac equation in Eq. (jBj), we have 

(i Hrfdn - m 7 c) ^ (x) = (259) 



or 



which yields 



since 



— 7 °-m 7 cjUS(*) = 0, (260) 



fko = rrijC , (261) 
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7. 5. Rotation of the wave functions 

The rotations of the wave functions follow from the discussion of Sec. 16.41 with an 
additional consideration of the vector k. On physical grounds, a rotation parameterized by 
the vector u of the state of a photon means rotation of the vector k into the vector k' , 
according to 



where R(u) is defined by Eq. fllOip . Similarly, the polarization vector is transformed by the 
spherical rotation operator 



k' 




(263) 




(264) 



and 




(265) 



where the rotation angle 6, axis direction u, k, and k 



are related by 



usm 9 



k x k'. 



(266) 



For applications, the rotation operator can be expressed as a function of k and k rather 
than u. From 




(267) 



one has for transverse polarization, A 



1,20 






1 + k' ■ k 



(269) 



and for longitudinal polarization, A = 0, 




(270) 



We thus have for a rotation of the wave functions characterized by the vector u 




(271) 



3 The identity r • k X k' = i(k' s fe| — fe s k'J) is useful here. 
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where k' ■ x = k ■ R 1 (it) x in the exponent of the wave function. This yields the result 



0. 



(272) 



according to the discussion in Sec. 16.41 

The expected transformation in Eq. (12711) can be confirmed by an explicit calculation 
based on the completeness of the wave functions. For a rotation of a transverse wave function, 
we write 



K(u) if>£l(Rr x (u) x) = I dx! S(x - x x ) K{u) ^[{Rr^u) x x ) e 



(k) / rt-l. 



— KlUlt 



E E / dfe i / ^;i 1 (^)<l| t 1 (^i)^W^|( J R- 1 H :Kl )e-^, (273) 



K\— >± Al=l 



where A = 1 or 2, and from Eq. (I164p and the subsequent remarks, it follows that the rotated 
wave function is also transverse, so Ai is restricted to 1 or 2. The evaluation requires the 
matrix element 



(«) / o-l, 



2(2vr) 3 

2 °KlK 



dxi 



Vl (fci) i? s (it) e A (fc) + e Al (&i) t • fci H s (ti) r ■ k e x (k) 



v — /tiifci-cei Kik-R 1 (u)xi 
At) t; 



4 (fc') i^H e A (fe) + e T Al (£') r-k'R s (u)r-k e x (k) S{kt - k') 



S KlK 5x 1 x8(ki - k'). 



(274) 



In the exponent in Eq. (12 74p . we have k ■ R~ 1 (u)xi = R(u)k ■ x± = k' ■ X\. The factor 
S KlK results from the requirement that k' — > k as the rotation angle 9 — > 0, i.e., that k does 
not change sign for an infinitesimal rotation. The factor 5\ 1 \ follows from Eq. (12641) and 
the discussion that follows it, together with the orthonormality of the polarization vectors. 
Substitution of Eq. (1273)1 into (12731) yields 



(275) 



in accord with the general arguments leading to Eq. (12711) . 
For a rotation of a longitudinal wave function, we have 



K{u) ^l{R-\u)x) = J dx 1 5(x-x 1 )TZ(u) ^ (R-\u) Xl ) 

= E [ dk ^ [ dx ^ ^(^l\^)n^^(R-\u) Xl ) (276) 
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where only longitudinal wave functions contribute, and the evaluation requires the matrix 
element 



d Xl ^l\ Xl )n{u)^{R-\u) Xl ) 
1 

= 5 KlK el{k')R s {u) eo(fe)5(fci-fc') 



dxi8 KlK ej( fc i) ^s(^) 6o(fc) e 



-/ciifei-asi /cifc-J? 1 (it)a;i 



= 5 KlK 5(k 1 - k'). 
Substitution of Eq. fl277j) into (12761) yields 

n{u)^ {R-\u)x) 
in agreement with Eq. 02711) . 



(277) 



(278) 



7. Velocity transformation of the wave functions 
7.6.1. Transverse wave functions 

Under the velocity transformation of a transverse photon by a velocity v, the four- vector 



k 



\k\ 
k, 



(279) 



is transformed to 

k' = V(v)k 



| k | ( cosh £ + £ ■ fc sinh £ 



\k\v r 



sinh £ + v ■ k (cosh £ — 1 ) 



and the wave function is expected to transform according to 



(280) 



(28i; 



The prime on the transformed function indicates that it is a function of modified polarization 
vectors, which in general are not simply rotated vectors corresponding to k — > k . The 
transformed transverse wave function is proportional to 



V(v) 



ex(k) ^ 
r ■ k i x (k) 



I + (r ■ 6) 2 (cosh( — 1) + r • vr ■ k sinh^j e\(k) 
t ■ v sinh + (r ■ v) 2 (cosh ( — 1)] r ■ k j e A (fc) 

r ■ k> e' A (fc') 



(282) 



where the fact that transformed wave function can be written in the form given at the 
right-hand end is based on the three identities: 

= cosh C + v ■ k sinh C, (283) 



I + (t • v) 2 (cosh £ — 1) + t ■ vt ■ k sinh £ \ e A (fc 
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[T ■ V 



(cosh £ — 1) + t ■ vt • k sinh £ j e\(k) 



(284) 



v sinh £ + [/ + (r ■ t>) 2 (cosh £ — 1)1 t • 

= r ■ fc' | J + (t • v) 2 (cosh£ — 1) + r • vt ■ k sinh^j 



(285) 



where 



k + v 



sinh ( + v ■ fc(cosh ( — 1) 



k' 

cosh C + v ■ k sinh ( 
Equation (I283j) determines the scalar multiplicative factor in Eq. (12811) to be 

£ = cosh ( + v ■ k sinh 



(286) 



(287) 



According to Eq. (12841) . the transformed polarization vector is in the plane perpendicular 
to k', that is, the transformed transverse polarization vector is also transverse. This is in 
contrast to the vector potential, which does not maintain transversality under a velocity 
transformation, so the Coulomb, or radiation, gauge condition is not preserved. The dif- 
ference is just the fact that the vector potential transforms as a vector, so that the angle 
between the space component of the vector potential and the space component of the wave 
vector is not necessarily preserved under a velocity transformation, whereas the polarization 
vector transforms as the electric field, i.e., as a component of a second-rank tensor, and 
Eq. (I284p shows that for this case the transversality is preserved. Equation (I285p shows 
that the lower components of the transformed wave function, as given in Eq. (12821) . can be 
written as r • k times the upper components in the same expression. This together with 
the relation k' ■ x = k ■ V~ 1 (v)x in the exponent of the wave function insures that the 
transformed wave function is a solution of the source-free Maxwell equation. It also follows 
from Eq. fl282l that if A 2 / Ai, then 



s't 

; A 2 



k') e Xl (k') = 0, 



(288) 



so the orthogonality of the polarization vectors is preserved by the velocity transformation. 

Equation (12811) can be be obtained by an explicit calculation based on the completeness 
of the wave functions, as for rotations. A velocity transformation of a transverse wave 
function is given by 

V(v)ipli(V-\ V )x) 

d Xl 8{x - xi) V(v) Tpl K {(V-\v) Xl ) 



«i->± Ai=0 



E / dfc i / ^(^^(^OV^)^^^ 1 ^)^), (289) 
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where A = 1 or 2, ti = t, and the primes on the wave functions indicate that the velocity 
transformed polarization vectors provide the basis vectors for Ai = 1, 2. For Ai = 0, 

d Xl ^(x 1 )V(v)^(V-\v)x 1 ) 
= ^ <W ej(fe') e' x (k') 5(k t - k') e- ik '° ct = 0, (290) 

= ^&-«Z r ■ k'e' x {k') 5{k x - k') e ik '° ct = 0. (291) 

In the exponent in Eqs. (12901) and fl29T|) . k-V- 1 {v)x 1 = V(v)k ■ %\ = k' ■ %\. As for rotations, 
the factors 5 +K and 5_ K result from the requirement that k' — > k as the velocity \v\ — > 0, 
and the last equalities follow from Eqs. (12841) and (12851) . For Ai = 1 or 2, we have 

xe -Klifci-£Ci e -Kifc-V r_1 («)xi 

= 6 KlK 6 Xl x£6(k 1 -k')e- Kik ' 0ct . (292) 
Substitution of Eq. fl292|) into fl289l) yields 

VW^ir 1 ^),) = (293) 

where 

fe^fc') = |/+(r-v) 2 (coshC-l) + T-i;r-fcsinhc}e A (fc), (294) 

in accord with Eq. (12811) . 

7. Longitudinal wave functions 

For a longitudinal solution, the four- vector /c in the invariant phase factor k ■ x, given by 

k = ( ° ) (295) 
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from Eqs. f)254l) - fl256p . has a zero time component. However, the transformed phase, given 
by k ■ V~ l (v) x = V(v)k ■ x = k' ■ x, where 



« = "W*=(*. + «f«*^c-i))" (296) 

does have time dependence. Thus the wave function is expected to transform according to 

EE*)^) = V{v)^ {V-\v)x) : (297) 

/c'->± A' = 

where the coefficients include the extra time dependence introduced by the velocity trans- 
formation. The sum over polarization states is necessary, because unlike the case of the 
transverse wave function, the transformed longitudinal wave function is mixture of longi- 
tudinal and transverse components. This is expected on physical grounds, since moving 
charges cause radiative atomic transitions. On the other hand, the transformed space-like 
wave vector does not the match the wave vector of either the longitudinal or transverse basis 
functions. Since the solutions are classified according to their three-wave-vector, there is a 
residual time dependence in the expansion that is included in the transformation coefficients. 

To be explicit, for A / 0, we consider specific polarization vectors. Let e[(k') be a 
linear polarization vector in the plane of k and v and e' 2 (k ) be a linear polarization vector 
perpendicular to k and v. These conditions, together with k'J e^(fe') = 0, yield (up to phase 
factors) 

e[(k') = , (298) 



1 - (y ■ k') 2 

iUk') T '* K , (299) 

'l-(v- k') 2 



where 



e' 2 (k') = r-k'e[(k'), (300) 
e'^fc') = T-k'e 2 (k'). (301) 

The transformed upper- and lower-component longitudinal wave functions follow from the 
expressions (for t = and x = 0) 

/ e (fe) \ ( [J+ (r • vf (cosh C - 1)] e (fe) \ _ ++( , ( e' (k>) 
V[V) { J ~ \ r-t)sinhCe (fc) 

i eiV(0) ( e[(k') \ g0)/ el(fe') . 

+ y/2 \T-k'i[Ck') ) + V2 { -t ■ k> e[(k>) ( - )U - ) 
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and 



K n io{k) J " " V [I + {t- v) 2 (coshC - 1)] e (fe) / ^ 00 W ^ e^fc') 



^2 V r-k'e' 2 (k') ) v/2 \ r ■ k' e' 2 (k') J ' 
based on the relations 

k' — v v ■ k' (1 — sech £) 



1 - (v- fc') 2 tanh 2 C 



k' s k'J [1+ (r •'0) 2 (coshC-l)] k = coshC\/l-(v-fc / ) 2 tanh 2 C ^ 
k' s k'^r ■ v sinh £ fc s = 0, 



(/ - k' s k$) [/ + (r • 6) 2 (cosh C - 1)] 4 



sinh ( tanh £ v • k'(v ■ k'k' s — v s 



I — k s k'j) r ■ v sinh ( k s 



1 - (v •fc') 2 tanh'C 

sinh £ r • fc'(i) • fc'fcg — v s ) 
1- (v ■ fc') 2 tanh 2 C 



with coefficients given by 



e o + o + (0)=e o T(0) = coshCVl-^-fcO'tan^C 



Terms that do not appear in Eq. (13021) or (13031) make no contribution 

£oV(o) = & + (o) = er _ (o) = ^iV(o) = & + (o) = e 2 _ o + (o) = o. 

An explicit calculation of Eq. (12971) is made by evaluating 

V{v)^l{V-\v)x)= j dx 1 5(x-x 1 )V(v)^(V- 1 (v)x 1 ) 



/ci-»± Ai=0 
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For Ai = 0, 



(2^3 / d*i ( <f(fci) ) V(„) ( *f } ) e-^-e-^W^ 
CoVW^^-feOe- 1 ^, (314) 



^ 1 ifcni ife'F 



=jiyj ix ^ ^Wu* le ""~' e 

= ^ 00 -(0)i(fc 1 -fe')e 1 *"' c ', (315) 



where 



fc' = + (316) 

^l + (ii-fc) 2 siiih 2 C 

and is the time component associated with the transformed space-like vector k 

k'° = vk sinh ( = v ■ k' tanh (. (317) 

For Ai = 1 or 2, 

x -Kiifcl-asi -ifc-V _1 (w)xi 

= 5 Kl+ 5 All £++(0) <5(fei - fc') e- ifc '° ct + 5 Kl _5 All er + (0) 5(fei + fe') e- ifc '° rf , (318) 



^g— Klifci-asigifc-y 1 ('") :E i 

= ^-^2^(0) 5(^1 - fe') e ifc '° ct + <Wai2&"(0) 5(^1 + fe') e ifc '° ct . (319) 
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Here both signs of Ki are included, because there is no continuity condition on the transverse 
solutions, which are absent in the limit of small velocity transformations. These results yield 
Eq. f)297p with the non-zero coefficients given by 

tio + (t) = e o + o + (0)e- ii> - fc ' tanhfrf , (320) 
= £oT(0)e ifi - fc ' tanhCct , (321) 
£++(t) = ^++(o)e i (l fe 'l-fi-fe'ta n hc) c t ) (322) 

£-+(*) = Q+( )e^l fc 'l +f) - fc ' tanh ^ rf , (323) 
&""(*) = ^(0)e i(l ^* tah()d , (324) 
&"(*) = g^e-^'l-^ - ^. (325) 

These transformed longitudinal solutions can be used, for example, to describe the fields of 
a moving charge by means of the expansion in Eq. (12481) . 

7. 7. Standing-wave parity eigenf unctions 

The parity operator changes the sign of the coordinates and includes multiplication 
by the matrix V = — 7°, so that the transformed wave function is also a solution of the 
Maxwell equations, as discussed in Sec. 16.61 We thus have 

W$0«0 = -TV®-*)- (326) 
With this definition, the parity operator commutes with the Hamiltonian in Eq. 



¥H = WP, (327) 
so we may identify eigenstates of both parity and energy. Since 

= <;!(*), (328) 

the parity eigenvalues are ±1. 

Transverse parity and energy eigenstates are 

{ ) ~ ^2^\Kir-ke x (k)smk-x )' ^ 



where 



v^t\^ = ±^t\x), (331) 

n^f\x) = Khe\k\^f\x), (332) 



41 



and A = 1, 2. These states are linear combinations of plane-wave states that form standing 
plane waves. Orthogonality relations are calculated with the aid of the integrals 

J dx cos k 2 ■ x cos fc x ■ x = Air 3 [S(k 2 — fei) + 5(k 2 + fei)] , (333) 

J dx sin k 2 ■ x sin fc x • x = 4n 3 [5(k 2 — fei) — 5{k 2 + fei)] , (334) 

J dx cos k 2 ■ x sin ki ■ x = 0, (335) 

which show that the states that differ only by the sign of k are not orthogonal. One of the 
overlapping states may be removed by including only states with k such that k ■ k > 0, 
where k is a fixed direction in space. With this restriction on k 2 and fei, the orthonormality 
of the states is given by 



dx^%f\x)^2C\ x ) = W^WCfc - (336) 

Completeness of these eigenfunctions, including the restriction on fe provided by a factor 
6{k • fc ), where the theta function is defined as 

{1 for x > 
\ forx = , (337) 
for x < 

follows from 

2 



E / dke(k.k Q )^j\x 2 )^M 

r~ V-L \ 1 " 

J dk0{k ■ fc )Z7" T (fe)(cos k ■ x 2 cosfe ■ x 1 + sinfe ■ x 2 sinfc • x±) 



K,TT^± A=l 

2 



(2tt) 3 

77 T (V 2 )5(a; 2 - : r 1 ). (338) 



Longitudinal parity and energy eigenstates are given by 

e (fe) sin k ■ x 


e (fe) cos k ■ x 




e (fe) cosfe ■ x 




^o' +) (^) = 


1 ( 




^-\*) = 


1 f 




^;o ,+) (*) = 


1 f 




^-\x) = 


1 f 


\/47r 3 \ 
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(339) 
(340) 
(341) 
(342) 



where 



Wjf^) = ^i^C*). ( 343 ) 
= 0. (344) 

As for the transverse eigenfunctions, there is overlap between states with opposite signs of 
k, so the same condition on k may be applied here, that is 9{k ■ ko) > 0, to eliminate the 
redundancy. With this condition on k 2 and ki, the orthonormality relation is 

J &x^^ ) \x)^^\x) = <WW(*a-*i)- 

(345) 

In addition, the longitudinal parity eigenfunctions are orthogonal to the transverse parity 
eigenfunctions. Completeness of the longitudinal parity eigenfunctions is given by 



E f dke(k-k )^\x 2 )^\ Xl ) 



K.TT— >± 



2 



(2tt) 3 

= n L (V 2 )5(x 2 - Xl ). (346) 
The combined completeness relation is thus 

E E / dfcflCfc-fcoJ^C^^W =Z*(* 2 -a: 1 ). (347) 

K,7T^± A=0 

7.8. Wave packets 

The plane wave solutions considered in this section are not normalizable as ordinary 
functions. Rather, integrals over products of solutions should be interpreted in the sense 
of distributions or generalized functions, like the delta function 26]. That is, they provide 
a well-defined value for an integral when they are included in the integrand together with 
a suitable weight, or test function. However, the plane waves can provide the basis for 
an expansion of a normalizable wave packet as a sum and integral over a complete set of 
solutions of the Maxwell equation. If f^ K \k) is a suitable function, we write 

9 f {x) = E/ dfc/^fc)^*), (348) 

K A 

and tff is a solution of the Maxwell equation 

-fdMx) = ^|dfe/| K) (fc) 7 ^^(,) = 0. (349) 
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Further, !?y will be normalized if is, because 



dx $j(x) &f(x) 



[ dk [ dk 'fi K> ( k ) I dx^\x)^ y ( X )fi?\k') 

= J2j dk \fi K) ( k )\ 2 = 1 - (350) 

K X 

In view of the orthonormality (in the generalized sense) of the plane-wave solutions, Eq. (1348P 
may be inverted to give 

fi K \k) = Jdx^(x)9 f (x), (351) 

where we have specified &f(x) = ^/(x)| t=0 , and the time dependence of the wave function 
is given by Eq. (13481) . 

An example is a normalized photon wave packet which at t — has (approximately) a 
wave vector k , a transverse polarization vector e 1 (fc ) ) and a Gaussian envelope of length 
a and width b centered at x = 0: 

9f(x) = , ) e^e-H/^ 62 ), (352) 

ah WJ \T-k e 1 (k ) J 

where xu = x ■ k k Q and Xj_ — x — x\\. For a and h large compared to |fc | _1 , the packet has 
a functional form that resembles a positive-energy transverse plane wave. From Eq. (13511) . 
we have 

W(U) = — ( — \ F {K) (k) e -[(fco-fcu) 2 a 2 /4+fci!>V4] 



/W(jfe) = _ J_ Ff ; (fc) e -^ K0 ~ KK| ^ ° / ^^° 'J , (353) 

2 \7T / 

where fey = fc ■ fco ^0; fej. = fe — fey, for A = 1, 2, 

F^(fc) = iei(fc)(/ + T-fcT-fc )6 1 (fco), (354) 

and 

F (+) (fc) = -Let^)^^), (355) 

Ft\k) = -^4(fc)T.feoei(fco). (356) 

The wave packet has small longitudinal components, because i*Q (fe) is not necessarily 

zero unless k = ±fe - It has negative-energy components, but they are also suppressed, 

particularly as a, b — > oo, because for k — > — , the exponential factors in Eq. (1353 j) favor k = 
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— k , and for A = 1,2, F^ K \— k ) = as compared to F^^fco) = Thus a cancellation 
between the upper-three and lower-three components of the wave function suppresses the 
contribution of negative-energy eigenstates to the wave packet. 
The expectation value of the Hamiltonian 7i, Eq. (I208p . is 

(V f \H\Vf) = -ihc J dx$}(x)a ■ V $ f {x) = hc\k \ = fiu , (357) 

and the expectation value of the momentum P, Eq. (I209p . is 

(Vf\V\V f ) = hk . (358) 
The initial probability density Q(x) is 

Q(x) = ^(-)^N = j,(^)%- 2 H/^i^), (359) 

with 

J dxQ(x) = 1, (360) 

and the initial energy density E(x) is 

E(x) = $\{x) U V f (x) = huj Q(x) - i hck ■ VQ(x). (361) 
1 2 

The real part of the energy density is proportional to the probability density for the photon, 
and the imaginary term, which vanishes upon integration to arrive at the expectation value 
in Eq. (135 7p . reflects the change in the initial probability density at the point x due to the 
motion of the wave packet. At a fixed point in the path of the wave packet, the probability 
density increases as the packet approaches and decreases after the maximum of the wave 
packet has passed by. The time-dependent probability density is 

Q(x) = &}{x)V f {x) = \e~ int/h W f {x)\ 2 , (362) 
and the change at t = is 

9Q{X) - £ »}(.)(« -**)»,(«> 

J Im V\(x)H & f (x) = -ck ■ VQ(aj). (363) 



dt 



t=o h 

2 



The gradient produces a vector that points toward the maximum of the wave packet, so 
that on the forward side of the packet k ■ VQ(a;) is negative and the probability density is 
increasing, as expected. Eq. (I363p also shows that the wave packet is initially moving with 
velocity c in the direction of ko- 
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The time dependence of the wave packet, Eq. (1352p at t = 0, is given exactly by Eq. (13481) . 
However, approximations may be made in order to obtain a more transparent expression. 
In view of the exponential factors in Eq. (J353]), the assumption that a, b 3> |fco| _1 implies 
k pa nk , and 

F^{k) pa F ( x K \Kk ) = S xl 5 K+ , (364) 
so that from Eq. (1535)) . !fy -> ^, where 

= ^flV/dfcf I 1 *. V-i(| fc M-^) e -[(^o- fe|l ) 2 aV4 +fe l6 2 /4]_ (365) 

This is a normalized positive energy wave function with polarization ei that is an exact so- 
lution of the Maxwell equation 7^(9^ $'f{x) = 0. Further simplifications are the replacements 
-+ fc in the polarization vector matrix and \k\ — > k ■ k in the exponent, which yield 
$' f -> with 

qr"(x) = ^ f 2 V ( \ /dfce- i ( fc ^ ct - fc - £C )e'[( fc0 " fc il) 2a2/4+fcife2/4 

1 8tt§ VW V r-k ei(fc ) / 7 

ei(fco) ^ \ e - i( ^- fc o-*) e -[(^MV+</& 2 ] (366) 
r ■ k ei(feo) / 

which is an approximate wave function with a normalized Gaussian probability distribution 

q"(x) = 2 / 2 y ^[(^^v^] (367) 

acr \7r / 

that moves with velocity c in the &o direction. 

7.9. Conservation of probability 

The formulation of the Poynting theorem in Sec. H] can be reinterpreted here to demon- 
strate conservation of probability. We define the probability density four-vector to be 

q »(x) = ~&(x)'fV{x), (368) 
where in the previous section Q(x) = q°(x). For the source-free case, S(x) = 0, Eq. (I5T)) is 

d^W(x)^^{x) = 0, (369) 
which is the statement of conservation of probability 

d 

— q°{x) + cV-q{x) = 0. (370) 
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Applying this relation to plane- wave states gives consistent, although trivial, results. We 
have 

= ^(x)^(x) = (2ir)- s , (371) 

which reflects the fact that the probability distribution for a plane wave is uniform through- 
out space and not normalizable. For transverse plane waves, A = 1,2, the probability flux 
vector is 

q(x) = ^(x)cx^ k %) = (2n)- 3 k, (372) 

and for longitudinal plane waves q(x) = 0. 
For the wave packet in Eq. f)352p . at t — 

q(x) = k q°(x), (373) 

so Eq. (I363p is essentially the conservation of probability equation evaluated at t = 0. 
Conservation of probability is valid for any solution of the homogeneous Maxwell equation 
from the definition of the probability density operator. However, it also happens to be valid 
for the wave packet represented by which is not an exact solution of the wave equation. 
The wave equation is not satisfied by because there is a small extra term resulting from 
the gradient operator acting on the perpendicular coordinate x±. On the other hand, q{x) 
for fy'j is proportional to k and ko ■ V x± = 0, so there is no corresponding extra term in 
V-q(x). 



8. Angular momentum eigenfunctions 

Radiation emitted in atomic transitions is characterized by its angular momentum and 
parity. In this section, wave functions that are eigenfunctions of energy, angular momentum, 
and parity are given; they are also classified according to whether they are transverse or 
longitudinal. The three-component angular-momentum matrices given here are to some 
extent parallels of the three- vector functions of (27| . 

8.1. Angular momentum 

The spatial angular-momentum operator is given by 

L = xxp = -ihxxV, (374) 

and following the example of the Dirac equation, the total angular momentum is given as a 
3x3 matrix by ^ 

J = L + hT, (375) 

where it is understood that the first term on the right side is a 3 x 3 matrix with L for 
diagonal elements and zeros for the restQ (In order to adhere to convention, we denote 



4 In some works, where electric and magnetic fields or the vector potential are three-vector valued fields, 
the spin operator is represented by a cross product. For example, Corben and Schwinger [27[ write J z & = 
L Z <S> + i e z x <&, where $ is a vector potential, and in Edmonds [28| , spin is represented symbolically as i e X . 
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both the current three-vector and the angular-momentum matrix by J. In either case, the 
meaning should be clear from the context.) The extension to a 6 x 6 matrix is 



J = ( q j^J =xxV + HS (376) 

and we have 

[H,J] = 0, (377) 

so eigenfunctions of both energy and angular momentum may be constructed. The vanishing 
of the commutator follows from the relations 

[t • V, L] = -ihr x V, (378) 
[r • V, r] = ir x V, (379) 
[r-V,J] = 0. (380) 

It is of interest to see that 3 commutes with 7i only for the (relative) combination of L and 
r given in Eq. (13751) . To obtain eigenstates of the square of the total angular momentum 
^J 2 and the third component of angular momentum J7 3 , given by 

J 2 = ( J q j 2 ) (381) 



and 



j* = ( ; 3 ) , (382) 

we construct matrix spherical harmonics that are analogous to conventional vector spherical 
harmonics and are three-component extensions of the Dirac two-component spin-angular- 
momentum eigenfunctions. Orthonormal basis matrices are given by 

e«= ( V 6<°>= 111, &-V= I I . (3831 



and they satisfy the eigenvalue equations 






r 2 



: 2&\ (384) 

where the 2 may be regarded as s(s + 1), with s = 1 as the spin eigenvalue. The matrix 
spherical harmonics are 

Y]P(x) = ^{lm-vlu\llj m)Y™~ v {x) &\ (386) 
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with vector addition coefficients and spherical harmonics in the notation of 28]. The spher- 
ical harmonics satisfy the eigenvalue equations 



L% m (x) = h%l + 1) Y^ix)., 



L^Y^x) 



hmY^ix) 



(387) 
(388) 



and by their construction, the matrix spherical harmonics are eigenf unctions of the total 
angular momentum: 



J 2 Yj?(x) 
J 3 YJP(x) 



h 2 j{j + l)Y™{x) 



HmYJP(x). 



Explicit expressions in terms of spherical harmonics are 

/ 



{j+m)(j+l-m) ym-1 ( - \ 

2j(i+i) r i W 

m Y m (x) 



{j-m)(j+l+m) V m+1 



2i0'+l) 



Y™+\x) 



(389) 
(390) 



(391) 



YJj +1 {x) 



( / (j+l-m)(j+2-m) 



Y, 



771—1 , 



(2j+2)(2j+3) "j'+l 



X 



(j+l-m)(j+l+m) 
(j+l)(2j+3) 

(j+2+m)(j+l+m) V m+1 



17?! (4) 



(2j+2)(2j+3) 



(392) 



/ / (j-l+m)(j+m) v m-l(~\ \ 
(2j-l)2j Y J-1 \ X > 



These functions are orthonormal 



(2j-i)j r j-A x ) 



V (2j-l)2j ^-1 W / 



which follows from the relations 



d«i^*(A)i^ i (*) = «y Wl «s mami) 



(393) 



(394) 

(395) 
(396) 



/^(^ 1 J2 m K m — v 1 z/ )(^ ra ~ v 1 v\l 1 ji m) = 5 
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(397) 



and they are complete 

J2YK l ^2)Y^(x 1 )=I5(cos9 2 -cos9 l )5( ( l> 2 -4> 1 ), (398) 

jlm 

based on the relations 

m — v 2 1 v 2 \l 1 j m)(l 1 j m\l m — v x 1 Vi) — 5 U2V1 , (399) 

3 

e M e Mt = /, (400) 

V 

Y,yr{x 2 )Yr\x x ) = 5( CO s9 2 - cos 00 ^ - 00, (401) 

where 6*j, 0j are the spherical coordinates of x^. 

An alternative set of matrix angular-momentum eigenf unctions is 

X{ m {x) = -^==L s Y 3 m (x), (402) 



Xt{x) = — 7 ^==r.xL s Yp(x), (403) 

Xt(x) = x B Y™(x). (404) 

For j = 0, X$°(x) = X$°(x) = 0. From a comparison of Eqs. fl39H -fl393l) to Eqs. (1402]) - 
d3D3D, one has 

X( m (x) = Y%(x), (405) 



= -J^^i^-J^r^-iC*). ( 406 ) 



*2 m (*) = -y^TT^i^ + y^TT^i^)- ( 40T ) 

In view of the relations in Eqs. (I405p - (l407p . the functions X^ m (x) are eigenfunctions of J 2 
and J 3 with 

J 2 Xr(x) = h 2 j(j + l)Xt(x), (408) 

J 3 X/ m (x) = hmXi m {x). (409) 

This can be confirmed directly from the definitions in Eqs. (14021) - (I404p with the aid of the 
commutation relations 

[L\U] = \he %jk L\ (410) 

[L\x j ] = ihe l]k x k (411) 
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and the tau matrix identities in Sec. El which provide the operator identities 

J l L s = L S L\ (412) 

[J\r-x] = 0, (413) 

fx s = x s L\ (414) 

These functions are orthonormal 

JdOX^(x)X^(x) = 5 i2h 5 j2jl 5 m2mi , (415) 

and they are complete 

J2x{ m (^2) Xi m \x x ) = I5(cos6 2 - cos 00 5(0 2 - 00, (416) 

ijm 

where the latter fact may be seen from the completeness of the matrix spherical harmonics 
and the relation 

j2xr( X2 )xr\ Xl ) = y, y ^2)y-\ Xi ). (4i?) 

i I 

The parity of the eigenfunctions is given by 

xn-x) = (-iyxr(x), (4i8) 

Xt(-x) = (-iy +1 Xn*), (419) 
Xt(-x) = (-ly^Xn*), (420) 
which follows from Y^ m (— x) = (— 1)- J l^ m (£). 

8.2. Helicity eigenstates 

Special cases of the transverse plane-wave solutions in Eq. f)217p are circularly polarized 
states with the polarization vectors in Eq. f)216p . They can be grouped with the longitudinal 
plane- wave states in Eqs. (12371) and (I238j) with the polarization vector in Eq. (12361) . These 
polarization vectors are summarized here as 

^i(e 3 ) 

where we have changed the label for A from 2 to —1 for this section. 

The states have a well-defined helicity; they are eigenfunctions of the operator for the 
projection of angular momentum in the direction of the wave vector • k p, |29j. In view 
of the relations 

L-ke ±ikx = 0, (421) 
r-ke x (k) = Xe x (k) (422) 

for the polarizations considered here, we have 

J-k^l(x) = Xh^l(x) (423) 

for these states. 
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i 


h e (e 3 )=| 


1 


; <U(e 3 












8.3. Transverse spherical photons 

Transverse spherical wave functions are given by 



V2 I -# T .V/ Wij (r)Xr(^ 



J_ / -T-Vf^(r)XrW 
V2 \ Kif^{r)X[ m {x) 



(424) 
(425) 



where r = |ac| and j > 1. They are transverse because Vj L s = Vj r • V = 0, so that 



, T(/t,7r) , 



n h {v)^,. ix) 



0. 



The wave functions are eigenfunctions of angular momentum with eigenvalues 
Eq. dSSOD] 



i,jm 



,T(K,7r), 
ui,jm 



h~mip u ,V.. (x) 



and they are eigenfunctions of H, with eigenvalue nhw 



T(k,7t) 



provided 



which is true if 



(426) 
(427) 

given by [see 

(428) 
(429) 

(430) 
(431) 



r <9r 2 



0. 



Solutions of Eq. (14321) are spherical Bessel functions given by [30 

jj(ur/c) 



fu;,j(r) oc 



hj(u>r/c) 



We employ the normalized solution 



to 



2 \C 



jj(ur/c) 



(432) 



(433) 



(434) 



for the wave functions; any other linear combination of spherical Bessel functions (with 
j > 1) is not integrable as r — > 0. The parity of the wave functions is 

V&fH*) = (435) 
V^5m _> (*) = (-l)^^^)- (436) 

This provides the conventional parity and angular-momentum attributes for electric and 
magnetic multipole radiation. Namely, ^uim( x ) * s magnetic 2j-pole or Mj radiation and 



i^jm \ x ) * s electric 2j-pole or Ej radiation. 

Alternative forms for the lower three components in Eq. (14240 or the upper three com- 
ponents in Eq. (14250 are obtained by writing (see [D]) 

r-V = l^rT-x+±-(L s xi + x s Li), (437) 

which yields 

r ■V/, 3 (r)Xr(£) = \ ^rf uJ (r)X> m (x) - ^ + 1} f u j{r)Xp{£). (438) 
Relations among spherical Bessel functions provide 

1 9 rf UJ (r) = ^[(j + l)/ u ,H(r)-j'UMi (439) 



r dr c 2j + 

— fw.jv) I i 

r c 2j + 1 



zAiW = TSTTT^-iW+W'')], (440) 



which together with Eqs. (14061) and (14070 yield the second alternative form 



°-t • V f u , 3 {r)Xt{x) = -JljLLf^r) Y^(x) + J ^A,, i+1 (r) (441) 

This latter form is useful in calculating the wave function orthonormality and completeness 
relations. 

An analogous longitudinal function is obtained by writing [Eq. (15721) ] 

d 1 

V s = x s — -— txL s (442) 
or fir 



and 



* Ft{x) = V s / UiJ (r)r(i) 



c 



£ L^^Xrm - VjU r +1) UAr)Xr(x), (443) 
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so that 



uo 1 



From the additional relation 

d 

together with Eqs. (I406p and (I407p . one has 



J 



3 + 1 



The orthonormality of the transverse wave functions is given by 



which takes into account the integral 



drr 2 f U2j (r) U uj {r) = <y(w 2 - wi). 



'•^W& w)t (*i) = n\v 2 )5{x 2 - Xl ) 



The completeness relation for the transverse wave functions is 

poo 

/ 

Jo 

K-rrjm 

which is shown in some detail by writing 

KIT 

where 

sr(x 2 , Xl ) = f^xr^f^xr^x^ 











and 



= i7j(V 2 ) [/. )i (r 2 )^(^)/ ( , )i (r 1 )^f(* 1 ) 
+/ w j-i(r 2 ) y3?_i(*2)/ w j-i(ri) 



+/ w j+i(r 2 ) lS+i(* a )/ w j + i(ri) ^(Ai) 
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which gives 

poo 1 

/ duY,Sr(x 2 , Xl ) = nJ(V 2 ) 5(r 2 - ri ) 

IT^Va)^-*!), (453) 



based on 



dw/ W j(r 2 )/ w j(ri) = 5(r 2 -ri). (454) 

r 2 r! 



Longitudinal spherical photons 
Longitudinal spherical wave functions are 







<••!::„» = r( ' I- (456) 



' ! V s <^(r)Y/*(a-) 
V s ^(r)^ m (x) 



It follows from the identity r ■ V V s = that they are longitudinal 

^ L (V)0*) = (457) 
n T (V)^(x) = (458) 

and that they are eigenfunctions of H, with eigenvalue 

WVffi(«) = 0, (459) 

with no condition on g. However, in order to have a complete set of longitudinal wave 
functions, a set of functions, indexed by the parameter k is specified here. The form of the 
plane-wave longitudinal solutions in Eqs. (1237P and (1238 j) and of the expansion of a plane 
wave in spherical waves suggest the choice 

g kd (r) = k^-j^kr) (460) 

for the radial wave function, where A; is a free parameter. This set of functions provides an 
infinite orthonormal set of degenerate (u = 0) basis functions for each j. 

From the form of V a in Eq. (1442 p and the fact that X 3 2 m (x) and X^ m (x) are eigenfunc- 
tions of angular momentum, it follows that the longitudinal spherical wave functions are 
also eigenfunctions of angular momentum 

J 2 & x ) = + (461) 

^V£S&(«) = *rniftM(x). (462) 
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They have parity given by 

= (463) 



«P<S(«) = (-lyVftmC*)- (464) 

With spherical Bessel functions for the radial wave functions, we have, following Eqs. (14421) 
to (gll, 



^0u(r)17(*) = W^j^ <7^-i(r) (465) 

This identity facilitates calculation of the orthonormality relation for the longitudinal wave 
functions, which is 



dx^kSLi^ku^mSv) = ^Tri^fWi^-fci), (466) 
with 

/•oo 

/ drr 2 g k2j (r) g kuj (r) = 5(k 2 -k 1 ). (467) 
Jo 

The completeness relation for the longitudinal wave functions is 

POO 

/ Mj2&*2)ti!$!M = n L (V 2 )5(x 2 - Xl ), (468) 

Jo _„•„. 



which follows from 



where 



T/"'(. E2 .x,) ^ i (r 2 )p feJ (r 1 )^(^ 2 )^(^ 1 ), (1701 



t 



and 



Thus from 



V 2 ff tj (r)r(£) = -A^iM^)" (471) 



dk gkj(r 2 )gkj{ri) = —Sfa-n), (472) 
o r 2 n 
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we have 



dkJ2 T t( x 2,x 1 ) = n^(V 2 )S{x 2 - Xl ), (473) 

jm 

which provides Eq. (j468]). That result, together with Eq. (1449 j) gives the full completeness 
relation. 

As an illustration of a role of the spherical functions, we revisit the example of a point 
charge at the origin, for which 

= ( v ;' ) • <™> 



In view of the integral 



roc ^ rzz 2 

jf d* s , w (r) = (475) 



one has 

/"OO 

r,L(-) 



V2 7re Jo 

which is the analog, for spherical solutions, of Eq. (12481) for plane wave solutions. 

9. Maxwell Green function 

A solution of the Maxwell equation for the electric and magnetic fields &(x) given a 
specified current source S(x), as they are related in Eq. (j52|) 



six] 



can be found with the aid of the 6x6 matrix Maxwell Green function T>u{x 2 — given 
by 

2 . 

A=0 ^ 
2 . 

~ dk <A(^)^A(X1) C(t! - * 2 ), 

A=0 ^ 

(477) 

where ip^x( x ) * s gi ven by Eq. (12171) . (12371) or (12381) . and Eq. (I25ip . In view of the relations 

7^ 2 ^g(x 2 ) = 0, (478) 

1 ^B{t 2 -t 1 ) = To^-cti), (479) 

1 ^6{h-t 2 ) = -70^2-^0 (480) 
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and the completeness of the wave functions, we have 

^d%V M (x 2 - Xl ) = X5{x 2 -x x ) (481) 

and 

V u {x 2 - x^-f^dl = -15(x 2 -x 1 ), (482) 

where 

5(x 2 — xi) = 5(ct 2 — cti) 5(x 2 — Xi). (483) 
In terms of the Green function, a solution for the electric and magnetic fields is 

$(x 2 ) = [ d 4 x 1 V M (x 2 -x 1 )S(x 1 ), (484) 



as is confirmed by the application of jfj,d 2 to both sides. In Eq. (14841) d 4 x\ = cdtidx\. A 
separation into transverse or longitudinal solutions may be made by restricting the sum over 
polarizations to A = 1, 2 for a transverse solution or A = for a longitudinal solution. 

The Maxwell Green function also may be written as an integral over the four-vector k 
of the plane-wave solutions. For this it is useful to make the separation into transverse and 
longitudinal components. For the transverse part we have 

A=l J 

- E f dk v&M^tw eMt2 ~' l} - 

A=l J 

and we employ the identities 



(485) 



n .,— iko(ct2— cti) 



e( t 2 _ h ) = J_ / dfc -g— — — , (486) 

kail + 16) - oj / c 



l f°° q— ifeo(c*2— cti) 



_ e Mt 2 - tl) 9{fi _ h) = dkQ (487) 



where the limit <5 — > + for the integral is understood. We also have 

a.k^ x (x 2 ) = \k\^(x 2 ) = ^$(x 2 ). (4? 
Together these relations yield [see Eq. (1228 j) 



2 „ 

P£(a*-*i) = ^E d * k 

71 A=l ^ 



J " r e -i fc o(ci2-cti) _. + . 
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p -ifeo(c*2— Ctl) _l\ 



fc (l + i5) + a ■ fc 

pifc-(sB 2 -asi) / / £,\2 q 

d 4 fc e o(c * 2 ~ c * l) 



(2tt) 4 7 7 °A;o(l + i(5) -7-fc V (r • fc) 2 

= W™/^^^ (489) 

where d 4 /c = dko dfc, and Cf indicates that the contour of integration over ko is the Feynman 
contour, which passes from — oo below the negative real axis, through 0, and above the 
positive real axis to +00; this is equivalent to including the factor (1 + 18) multiplying fc in 
the denominator and integrating along the real axis. 

For applications, it is useful to consider an alternative form for the transverse Green 
function. Taking into account the relation 



n T (k) - - n T (k) ^° fco ~ 7 ' fc 

1 V*;o(l + i5)-7-fc ~ } k 2 -k 2 + i8 

k (r • k) 2 —t ■ k 



T k ~k (T- k) 2 J k 2 + 18 

we have, from Eq. (J489]), 

m{2 l} - (2tt)4 dko 6 ^ -ir-V -k nJ(V) J J dk k 2 

dk e 

dk e - iko(ct2 - ctl) 



(490) 



k 2 + iS 



2 /"OO 



^■f) 2 -(fcg + i^/ 2 r.f\ eW /2 H 
(fc 2 + i5) 1 / 2 r-f -fc (r-f) 2 J \r\ ' 1 j 

where r = x 2 — aci, the gradient V is with respect to r, and the branch of the square 
root in the exponent is determined by the condition Im(fcg + i<5) 1//2 > 0, which specifies that 
(&0 + i5) 1//2 — ^ I A^o I for real values of ko. In the last line of Eq. (14911) . higher-order terms in 
(ko M) -1 are not included, but the exact expression follows from the formulas in [El 
For the longitudinal Green function, we write 

V^(x 2 - Xl ) = Jdk^ \x 2 )^ \x 1 )e^- c ^6(t 2 -t 1 ) 

~ /dfc V$(*a$$(*i) e £(ct2 " cil) - fc), (492) 

where damping factors with e > are added so that the Green function falls off for large 

time differences. In addition, we assume that the longitudinal wave functions are solutions 
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of the Maxwell equation with an infinitesimal mass m e included, as given in Eq. (12591) . in 
order to be able to use the Feynman contour to specify the path of integration over ko in 
relation to the poles of the integrand. We employ the identities 



g-irrijC 2 '{t2-t\) i '% e ~e(ct2~cti) 



_ Am t c 2 (t2-t\)/h e(ct2— cti) 



0(h 



i 

2^ 
i 

2^ 



e -ifco(ct2-cti) 

dfc 

-oo k + le 

oo 

dfc , ,. 
k — le + m € c/n 



-ifc (ct 2 -cti) 



and 



7°<J(* a ) 



o 



(493) 
(494) 



(495) 
(496) 



to obtain 



e -ifco(c*2-cti) 

fco + ie — m e c/h 

p -ifeo(ct2— cti) 



(+)/ 



-t(«) 



1 
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d 4 £; 



A; — ie + m e c/h 

g— ifco(ct2— cti) 



C F 



A; — y°m e c/h 



5^ ^2(^2)^0(^1) 



1 



(27T) 



i7 L (V 5 



3 -ifc-(x 2 -xi) 



d4fc — 1 TiT 

c P y^-meC/h 



(497) 



Here the limit e — > would be undefined without the mass term. A concise alternative 
expression for the longitudinal Green function is obtained by the substitution of the partial 

and /r ° 



completeness relations that follow from Eqs. 



into Eq. (TTO : 



V\ A {x 2 - xi) 



n L (v 2 )5(x 2 - Xl ) 



ie(h-K 








(498) 



The transverse and longitudinal Green functions in Eqs. (14891) and (149 7p differ only by 
the type of projection operator and the infinitesimal mass term in Eq. (1497p . However, such 
a mass term in the last line of Eqs. (I489p would not change the relation of the path of 
integration over k to the location of the poles of the integrand, so it could also be included 
in that expression. In particular, the poles in Eq. (14901) at ko — ± (k 2 — W) 1 ^ 2 would move to 
k = ± [fc 2 + (m e c/h) 2 — i<5] 1//2 . These poles lie on curves in the second and fourth quadrants 
of the complex k plane, whereas the Feynman contour passes through the first and third 



quadrants. Thus, we may write V^x-2 — x\) + T>\j i {x 2 

T> M (x 2 - x x ) 



(2vr) 4 . 
60 



d 4 fc , 



Xx) = T> M (x 2 

e -ik-(x 2 -x 1 ) 



Xx), with 



m e c/H 



(499) 



This result is a covariant Green function for the Maxwell equation which is of the same form 
as the well-known Green function for the Dirac equation. A formal coordinate-representation 
is 



V M (x 2 -xi) = — ^S(x2 -xx). (500) 

The fields given by Eq. f)484p represent a particular solution of the Maxwell equation. 
Any solution of Eq. (I484p for S(xi) = 0, such as the field of a static charge distribution, may 
be added to the particular solution, and the sum will be a solution with the same source 
function. In fact, even if the three- vector current density vanishes in the distant past and 
future, there could be a static charge distribution that persists, with or without a net total 
charge, for which the fields would be non-zero indefinitely. To deal with this case, we obtain 
an expression that takes into account the possible fields in the past and future by writing 
the time derivative 

dx 1 V M (x 2 - x^o^ixi) = / dx 1 V M (x 2 - £1)70 d 1 &(xx) 



d(ch 

+ J d Xl V u {x 2 - xi) 7o 9 1 ^(xi), (501) 

and for fields that vanish for large space-like distances, we write 

J dxi Vi ■ V M (x 2 - £1)7 $(xx) = J dxiV M (x 2 - £1)7 ■ V x ${x x ) 

+ J d Xl V u {x 2 - xi) 7 • Vi &{xi) = 0, (502) 

where the result is zero, because it may be written as an integral over the bounding surface, 
by the Gauss- Ostrogradsky theorem. The sum of Eqs. (15011) and f!502j) is 

d f 

g ^ ct j J dx 1 V M (x 2 - £1)70 

= J dx 1 V M (x 2 - x 1 )'j fl 'd^(x 1 ) + J dx 1 V M (x 2 - x^^d^ixi) 

= - dx 1 5(x 2 - x 1 )^(x 1 ) + / dx 1 V M (x 2 - xt) E(xi). (503) 



Integration of Eq. (15031) over ti from ti to tf, where ti < t 2 < if, yields 

${x 2 ) = J dx 1 V M (x 2 - a;i)7o<P r (xi) ^ - J dx l V M (x 2 - x 1 )^(x 1 ) 



ti=t t 



+c£ dtx J dx 1 V M (x 2 - xx) S(xx) (504) 
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or 



V(x 2 ) = 




As a consistency check of this expression, we note that it properly reduces to the expected 
result for the field of a constant charge distribution. In this case, the current source term 
vanishes, the initial and final fields are the same, and they are purely longitudinal. As a 
result, only longitudinal functions with no time dependence will contribute to the sum over 
states, which is just the longitudinal completeness relation, and Eq. (1505 j) reduces to the 
proper identity. 



10. Applications of the Maxwell Green function 

The Maxwell Green function is used here to calculate the radiation fields of a point 
dipole source as an example of an application. Only the large distance transverse fields are 
considered, and they are given by 

^d(x 2 ) = J d i x 1 Vl l (x 2 - x-i) E d (xt), (506) 

with the source term 

5.M = (-"»^<;> e "" d '). (507) 

The classical source for dipole radiation is a charge q with position 

x d (t) = x cosu d t (508) 

which produces a current density 

j c \(x) = qS(x -x d (t) )x d (t) 

« —uj d d5(x)smuj d t, (509) 

where d = qx . This is the real part of 

j(x) = -iuj d d6(x)e- iuJdt , (510) 



which is the source current for the radiation [17j . For the transverse Maxwell Green function, 
we use the expression on the last line of Eq. (149ip . Integration over t\ yields a factor 
2ii5(k c — u d ), and evaluation of the integration over k follows. The result is 

c J dti T>m(x 2 - xi) E d {xi) 

fiock ( (r • r) 2 -r • f \( j s ( Xl ) \ e^ e _^ 2 + _ ^ (5n) 



4ni \ t f — (t ■ r) 2 J \ 
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where k = uj^/c. Since the source is point-like at the origin x± = 0, r — x 2 , and 

Ux) = *1 ( )^ e — + .... (512) 

47Te V T ' Xd * J \ X \ 

The time- average differential radiated power, based on Eqs. (|54|) to (159]) with a factor 1/2 
from the time averaging 1?J is 



^ = l -x 2 x-S(x) = C -^x*y £ (x) 1 -x>V £ (x) 

- S^^-^wk^-^-^ (513) 

which is the well-known result. 

A more realistic example is the radiation produced by a Dirac transition current, which 
is given by 

5b(«) = ( -"° c f <*> ) = ( *.*,(*) r*>> , f (514) 

where fa and 0/ are the initial and final hydrogen atom Dirac wave functions, here a is the 
4x4 Dirac matrix, and 

u if = — ^ (515) 



is the frequency corresponding to the energy difference of the transition. The factor of 2 
multiplying the matrix element accounts for the difference between a classical dipole moment 
and the quantum mechanical dipole moment operator in Eq. f)520p . (See the footnote on p. 
407 of El.) We have 



dt 1 Vl l (x 2 - xi) Sd{xi) 



ik ( {r-ff — t - v \ ( Wixx) \ e ife H 



47re c V r r -(r • ff J V 



where k = Uif/c. For distances far from the source atom, \x%\ ^> \x\\, k « r ~ £2, and in 
the exponent fc|r*| = /c|cc2| — k • X\ + . . . , which yields 

^d(^2) = J d 4 xi 2?m(2j -aci) ^d(^i) 

- JL / da, ( {T • ^ e- ifc - ^ e"^ + f 517) 

" A,e cJ dXl { r-fc&M J" \x 2 \ e +••• ' (517) 
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The average radiated power is 



d/p_ 



x\x 2 ■ S(x 2 ) 



^-x\ ^ D (x 2 )7 • kV D (x 2 ) 



= jdx 1 j^(x 1 )^(r-k) 2 Jd^jTix'Je-** 

= huj lS ^Y, [ dx4(x)e x (k)- a e ik - x Mx) ! daj'4(aj')6 A (fe).ae- w i (x'), 

^ A=l ^ ^ 

(518) 

where a = e 2 /47reo^c is the fine-structure constant. The radiated power integrated over 
directions of the vector k may be interpreted as fkVifAif, where Aif is the radiative transition 
rate for i — > f, that is, the probability that the atom providing the source current makes a 
transition from state % to state / in one second. This gives 



At 



if 



akc 



r 2 

/ dQ k (i \ex(k) -ae ikx /) (/ e A (fc) • a e 

J A=l 



-\kx 



(519) 



which is the same as the relativistic radiative transition rate given by QED (see [F]). In 
the dipole approximation e lkx —>■ 1, (i \ot\ f) — ik(i \x\ f), which follows from the identity 
[H, x] = [col ■ p, x] = —ihcat, where H is the Dirac Hamiltonian, and integration over k 
yields the familiar result 



.4 



Aau, 



if 



if 



3c 2 



</l*l<>| : 



(520) 



11. Summary 



In Eq. (1521) . two of the Maxwell equations, Eqs. (j2J) and ([3]), are written in the form of 
the Dirac equation without a mass, but with the addition of a source term S{x): 



where the gamma matrices are 6x6 versions of the Dirac gamma matrices in Eq. (j4"6"I) . and 



E s (x) 
icBJx) 



Z{X 



-fi cJ s (x) 




from Eqs. (|4"T1) and fl5Tl) . The source-free version of this equation, with S(x) = 0, can be 
written as a Schrodinger-like equation from Eq. (|45|) or 
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(521) 



where the Hamiltonian, Eq. (1208p . 

H = — i he a ■ V 

is the analog of the Dirac Hamiltonian for the electron. The factors of h are not essential 
here, but they are introduced to provide the conventional units of frequency and energy. As 
with the Dirac wave functions, where all four components are necessary to describe an elec- 
tron bound in an atom relativistically, all six of the components of the photon wave function 
apparently are necessary to properly account for the space-time properties of electromag- 
netic fields. There are three polarization degrees of freedom, two for radiation and one for 
electrostatic interactions, and relativistic covariance requires twice that many components. 
Alternatively stated, six complex functions are necessary to describe the six components of 
the electric and magnetic fields, and they are coupled by the Maxwell equation and Lorentz 
transformations. 

According to Eq. (15211) . as in Eq. (12581) . the time dependence of the solution is given by 

$( x ) = e - mt/h &(x). (522) 

The time-independent solutions may be expanded in eigenf unctions of the Hamiltonian with 
eigenvalues E n given by 

W n (x) = E n V n (x), (523) 

where n is a set of parameters that characterize the state represented by the wave function. 
For each eigenf unction, one has 

y n (x) = e - iEnt/h V n (x). (524) 
The eigenfunctions are orthonormal 

J da; $1 (x) (x) = S n2ni , (525) 

and they are complete 

£ ¥ n {x 2 ) Vl{x x ) = 8{x 2 - x x ). (526) 

n 

The state index n includes continuous variables, so Eq. (I525p has delta functions in those 
variables on the right-hand side, and the summation symbol in Eq. (I526p includes integration 
over those variables. 

States considered in detail in this paper are propagating plane waves in Sees. 17.11 and 
17.21 standing plane waves in Sec. IT. Tj. and angular-momentum eigenstates in Sees. 18.31 and 
18.41 The propagating plane-wave states are eigenfunctions of the momentum operator in 
Eq. (JUD 



V = -ihlV. 
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This operator commutes with the Hamiltonian, [TC, T > \ = 0, and eigenstates of both energy 
and momentum are given in Eqs. (12171) . (12371) . and (12381) . These plane- wave states are further 
characterized by polarization vectors given in Eq. (I210p or (I235p . Linear combinations of 
the traveling plane waves, combined to give standing-wave parity eigenfunctions, are in 
Eqs. (E2HD, ([MOD, and fl3^|) to (fM2|) . The angular-momentum operator, Eq. (13761) . is 

J = xxV + hS, 

where the spin matrix, Eq. (11541) . is 

T 
T 

and r is given by Eqs. ffl~5]) to ffTT|) . One has [H,&] = in Eq. (13771) . and simultaneous 
eigenfunctions of energy, angular momentum squared ^T 2 , third component of angular mo- 
mentum J7" 3 , and parity are given in Eq (14241) . (14251) . (14551) . and (14561) . All three sets of 
eigenfunctions listed above are shown to be orthogonal and complete, as in Eqs. (15251) and 



The eigenfunctions considered here are not normalizable wave functions. However, they 
provide basis functions for the expansion of a normalizable wave packet, as discussed in 
Sec. 17.81 For the sum 

V f (x) = Y,fn# n (x), (527) 

n 

from the orthonormality of the eigenfunctions one has 

f n = J dx*l(x)V f (x) (528) 

and 

/ dx $){x) V f (x) = \fn\ 2 = 1 (529) 

J n 

for suitably chosen f n . For the example of a Gaussian wave packet in Eq. (I352p . the expec- 
tation value of the Hamiltonian, in Eq. (I357p . is 

(V f \H\Vf) = hu , 

where ujq = c\k \ is the frequency corresponding to the wave vector k of the wave packet. It 
is clear that Eq. (I357p applies in more generality than just to the wave packet in Eq. (13521) . If 
the Gaussian shape function were replaced by any normalized real function, the expectation 
value of the Hamiltonian would still be exactly ftujQ. For the wave packet in Eq. (13521) . the 
expectation expectation value of the momentum operator, Eq. (13581) . is 



(& f \V\Vf) = hk , 
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and the expectation value of the projection of the angular momentum in the direction of the 
wave vector is 



Jk c 



!P>) = ne t 1 (fc )r-fc ei(fco), (530) 



where the result depends on the polarization state represented by ii(k ). For circular 
polarization, Eq. fl 2 1 6 [) . the expectation value is ±h, while for linear polarization, Eq. (12151) . 
it is 0. The real part of the energy density is Hujq times the probability density 

Re &}{x)H9f{x) = hu V}(x)V f {x) (531) 

for the wave packet. 

For any wave packet, represented by \Pf, the photon probability density four- vector 
defined in Eq. (I368p is 

q»( x ) = l/ f (x)^V f (x), (532) 
and the differential conservation of probability is given by Eq. (13701) 

^q° f (x) + cV-q f (x) = 0. (533) 

This is valid for any solution of the homogeneous Maxwell equation. By integrating over 
a closed volume and converting the divergence to an integral of the normal component of 
the vector over the surface, one obtains a statement of conservation of probability for the 
volume. If the surface of the volume is taken to infinity in all directions, where the wave 
function vanishes, this expression shows the time independence of the normalization of the 
wave function 







dx &Ux) & f (x) = 0. (534) 



Of course, this does not give meaningful results for a plane wave, because in this case, the 
probability density is constant over space and is not normalizable. 

For electromagnetic fields and photons, Lorentz invariance is a necessary consideration. 
In Sees. 16.41 and 16.51 it is shown that 



a ix 



where the primes indicate that the field and source have been transformed by either a rotation 
or a velocity boost. For rotations represented by the vector u = 9u, the transformed field, 
in Eq. ffT44j) . is 

V\x) = TZiuj^R^iujx), 

where R(u) is the coordinate rotation operator in Eq. fll(J3j) and 

K{u) = e-' lSu 
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in Eq. (I153p . The source term transforms in the same way. For velocity transformations, 
corresponding to the velocity v = ctanh^i), Eq. (11651) is 



v(v)y{y-\v)x), 



where V(v) is the coordinate velocity transformation operator in Eq. H 123[) and 



V(v) 



e CK-v 



in Eq. (I169p . where 



T 
T 



in Eq. (I170p . The transformation of the source term under a velocity boost is noteworthy. 
In the presence of a non-zero source, the Maxwell equation is invariant, but the left- and 
right-hand sides do not transform separately. As shown in Sec. I6.5[ the derivatives acting 
on the fields produce terms that combine with the original source term in such a way as to 
produce the velocity transformed source term, even though it is the three- vector current. In 
the absence of sources, the transformation reduces to a more conventional form. 

In the case of photon wave functions, sources are taken to be absent and the wave 
functions are solutions of the homogeneous Maxwell equation. The Lorentz transformations 
of the plane-wave functions are explicitly shown in Sees. 17.51 and 17.61 As with the Dirac 
equation for an electron, the eigenvalues in Eq. (I523p may be either positive or negative, 
and here they also may be zero, for both plane-wave and spherical-wave eigenfunctions. The 
negative eigenvalues, which are associated with relativistic invariance, are necessary in order 
to have a complete set of solutions satisfying Eq. (I526p . It is relevant to note that for the 
six-component wave packet in Eq. (I352p . there is an interaction between the upper-three 
components and lower-three components, evident in Eq. (I354p . that suppresses the role of 
the negative energy states. 

In Sec. El Eqs. fl86l) and (IBTl) . orthogonal transverse and longitudinal projection operators 
are defined: 



These operators commute with the Hamiltonian, the momentum operator, and the angular- 
momentum operator 




where from Eqs 



(HDD and flZU) 




[n, n} = [v,n} = [j, 77] = o, 
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(535) 



where IJ represents either projection operator, i7 T (V) or 77 L (V), so all the eigenstates 
considered in this paper are classified as being either transverse or longitudinal, with 

n T (V)^(x) = ^(x), (536) 
i7 L (V)^ L (x) = ^ L (x), (537) 

respectively. The transverse states describe radiation and have non-zero eigenvalues in 
Eq. f)523p . while the longitudinal states correspond to electrostatic interactions, with an 
eigenvalue of zero. An exception is that the longitudinal states may have a non-zero eigen- 
value if a hypothetical mass term is considered, as discussed in Sec. I7.4L The projection 
operators commute with rotations, but not, in general, with velocity boosts. However, as 
shown in Sec. 17.6.11 the velocity transformed transverse plane-wave states are also transverse. 
This corresponds to the fact that radiation may be treated relativistically independent of 
electrostatic interactions. On the other hand, as shown in Sec. l7.6.2l the velocity transformed 
longitudinal states have both longitudinal and transverse components, corresponding to the 
fact that moving charges may excite radiative transitions. Both transverse and longitudinal 
states are necessary in order to have a complete set, as in Eq. (15261) . 

A solution of the inhomogeneous Maxwell equation may be obtained with the Maxwell 
Green function, as discussed in Sec. [9j The Green function satisfies the equation 

Jpd% V u {x2 -xt) = 15(x 2 -xi) 

in Eq. (j481|) . and a solution of the Maxwell equation is given by 

&{x 2 ) = J d 4 XiP M (^2 - xt) E(xt) 

in Eq. (14841) . It is shown, by summing over the complete set of plane- wave solutions, that 
the Green function is 

j r Q—ik-(x2—xi) 

V M (x 2 - xt) = — — / d 4 /c — — 

(2tt) 4 J Cf - m e c/h 

in Eq. (J499]), which is the same form as the Dirac Green function, except that here it is a 
6x6 matrix instead of a 4 x 4 matrix. In this equation, Cf is the Feynman contour and 
the infinitesimal mass is included to resolve an ambiguity in the longitudinal contribution. 

In SeadOJ applications of the Maxwell Green function are made, including a calculation of 
radiation from a Dirac-electron current source. In this example the six-component Maxwell 
formalism couples radiation to the Dirac current relativistically with a result that is the 
same as the result of a calculation that starts from Feynman-gauge QED. 

12. Conclusion 

We conclude that the criteria for properties of a single-photon wave function proposed 

in the introduction are met by the formalism described in the subsequent sections. In 
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particular, the example of a photon wave packet provides a normalizable solution of the wave 
equation whose properties can be verified by explicit calculations. It yields the unanticipated 
result that for virtually any probability distribution, under rather mild assumptions about 
the form of the wave packet, the expectation value of the Hamiltonian is exactly ( H ) = hw , 
where u> is the frequency associated with the wave vector of the packet. 

A. Velocity transformation of electromagnetic fields 

The velocity transformation of electromagnetic fields is derived here without invoking 
potentials for completeness. With the aid of the identity Vjr ■ cB = (V x cB)J , Eqs. (JTJ 
and fl2j) may be written as 

VjE c = p c 2 p, (538) 

<9_E7 T 

-^-V T c f-cB = -fi cJj (539) 

or 

dJgF = fi J T , (540) 

where 



1/0 -E, 



T 
c 

r \ E r f cB 

is the field tensor [see Eq. (13*81 ] and 



(541) 



J = ( °P ) . (542) 



J, 

Since V(v) g V(v) = g, Eq. (15401) is equivalent to 

dJV(v)gV(v)F(x)V(v) = fi J T (x)V(v) (543) 

or 

dJgV(v)F (V- 1 (v)x)V(v) = ^J T (V~\v)x)V{v). (544) 
Assuming the current transforms four- vector, 

J'( x ) = V(v) J{V~ l {v)x) , (545) 
Eq. (I540p will be invariant if the field tensor transforms according to 

F'(x) = V{v)F{V-\v)x) V(v). (546) 
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Direct calculation yield^f] 



F'(x') 



1 

c 



E' c {x) 



f • cB'(x) 



-E'J(x) 



where x' = V(v) x and 

E' = E c cosh ( — v c v ■ E (cosh ( — 1) — f • vcB c smh(, 
cB c cosh Q — v c v ■ cB (cosh ( — 1) + f ■ v E c sinh (. 



c 

cB' 



(547) 



(548) 
(549) 



These relations are equivalent (up to the velocity sign convention) to the electric and mag- 
netic field transformations in 17j, and in the spherical basis they are 



E' 



E s + (r • v) 2 E s (cosh( — 1) +ir • vcB s sinh(, 
cB' s = cB s + (r • v) 2 cB s (cosh( - 1) - ir • vE s smh(. 

The Cartesian transformation in Eqs. (15481) and (I549p can be written as 



where 



and 



K 



%{x) 



f 
-f 



E c (x) 

cB r (x) 



X 



(550) 
(551) 



(552) 



(553) 



(554) 



Similarly, the Cartesian form of the rotation transformation, corresponding to the operator 
in Eq. (jT55]l . is 







= e s - u V c (x), 


(555) 


where 










s = 


my 


(556) 


The components of the matrices in 


Eqs. (1553p and Q556j) have the commutation relations 






- c S k 


(557) 






£~k 


(558) 






qk 


(559) 


characteristic of the Lie algebra of Lorentz transformations. 




5 The identity e ljk = ei jk v l v l + e ak v j v l 


+ e i3 iv k 


v may be useful here. 
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B. Inverse Laplacian 

For cases where the integral definition of the inverse Laplacian converges poorly, we use a 
generalized definition that includes a damping factor to resolve ambiguity in the intermediate 
steps of the calculation. From the equation 

(V 2 -e 2 ) 6 = -4ir8(x-x'), (560) 



\x — X' 



one has 



1 1 p -e[sc-a>'| 

s ( x ~ x ') = -jzjzr-jy- ( 561 ) 



V 2 - e 2 An \x - x' 

Multiplication by f(x') and integration over x' yields 

1 1 r p-e|a!-as'| 

— -f(x) = -- dx f - -f(x'). (562) 

V 2 — e 2 An \x — x'\ 



We thus have, for example, 



by direct calculation of the integral. 



C. Coulomb matrix element 



The calculation of the Coulomb field matrix element in Eq. (12481) requires evaluation of 
the integral 

dx ^ e -^ = lim /dae-eM^e-** 



|£C| 3 J | a; | 3 

An 

- (564) 

which is defined with the convergence factor e" e|3:| . The inverse transformation requires the 
integral 

dk^je ik x = lim [ dke~^4ie ik - x 
\k\ e^O J \k\ 

\xr 



which confirms the result that 

[ 1 [ai.1j,(+)^ - ? ( x * 



/ d*j^|V$(aO 



^ftj \k\ rfe ' uv ' 47re |aj|3 V 

(566) 
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D. Separation of the transverse and longitudinal gradient operators 

The transverse gradient operator r ■ V is separated into radial and angular parts by 
writing 

T-V = T ■ XT ■ XT • V + X s x\t • V 

q I 

= -^r-x-T-xV s xl + —x s Ll (567) 

where the line over the gradient operator indicates that it does not act on the unit vector 
directly to the right. That term is 

t ■ x'V.xI = t ■ x I V s ^l + - x s xt — - I 



Thus 



* > — ' \ v S^s 1 •"S'^s 



— — L s x\ t ■ x. (56£ 

nr r 



Id 1 

t-V = - — rr- x + — (L s xl + x s Ll) . (569) 
r or nr v ' 



Acting on L s , the transverse gradient operator yields 

Id _ _ I 



t-VL s = -frr • xL s - -^x s L 2 , (570) 



where I/J£ 8 = — L 2 . 



For the longitudinal gradient operator V s , the identity 



provides 



d i 

V = X—-—XXL (571) 
or nr 



d 1 

V s = x s — -— txL s . (572) 
or nr 



E. Exact transverse Green function 

The exact transverse Maxwell Green function follows from Eq. (149 ip together with the 
relations 

e -w\r\ e -HH / \ \ 

t-V—— = -WT-r—— (l + — — (573) 



w\r\ 



and 



i _— to|t*| i 

^ V = "S / <iX T^°— = "db (1 - <=-"'-') • < 574 > 



1 


e 


-w\ 


x\ J 




r — a? 
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x\ 


If 2 


r 



which yield 



V l V J ' e~™ |r| 



r yi e -uj|r| 



+ 3 



-w\r\ 



-w\r\ 



w\r\ 



+ 



u> 2 |r| 3 u> 2 |r| 3 



and 



-iu|r| 



t ■ r 



+ [(r-f) 2 -2f s f s t ] 



e -w|i"| g— iu|r| 
^^7 + 



u> 2 |r| 3 u> 2 |r| 3 



(575) 



. (576) 



F. Radiative decay in quantum electrodynamics 

In QED, the radiative decay rate of an excited state may be obtained from the imaginary 
part of the radiative correction to the energy level of that state 



f 



-2Im(A£;) 



(577) 



where the sum is over all states with a lower unperturbed energy. This gives a correction 
to the level that, roughly speaking, results in an exponentially damped time dependence for 
the population of the state: 

| e -iAEt/R|2 = e - (578) 

For one-photon decays, the rate is included in the second-order self-energy correction to the 
level. An expression derived from Feynman-gauge QED that includes some of the real part 



and all of the imaginary part of this level shift in hydrogen-like atoms is given by [31 



ah 2 c 2 

1~ 2 Jfu L 

ah 2 c 2 



AE, 



dk 



1 



8, 



.a 



k j k l 



a j e ik ' x - 



4tt 2 



Jhck<Ei K A=1 j N 1 



H - Ei + hck- i8 

kx 



a l e~ ikx 



f 



x- 



/ e A (fc) • ae" 



ikx 



(579) 



E f -Ei + hck - iS 

where H is the Dirac Hamiltonian. The integrand is real, except for the imaginary infinites- 
imal in the denominator, for which 

1 



Im 



E f - Ei + hck - i5 



ir5(E f -Ei + hck) 



and hence 



E^/ = E ? / d ^E | • oeife '1 (/ | e ^^) ■ ae 

/ / 11 J A=l 



ikx 



(580) 



(581) 



with the restriction < Ef < Ei on the sum over /. The contribution to the decay rate 
from each final state / coincides with Eq. f)519p for the transition rate Aif. 
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